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ABSTRACT 


This  report  describes  an  analytical  model  for  predicting  the  emission  of 
radiation  from  a  jet  plume  in  the  mid-infrared  spectral  region.  It  is  assumed 
that  the  dominant  radiation  arises  for  the  CO2  molecule.  Results  are  therefore 
reported  for  the  4. 3-micrometer  band  of  gaseous  carbon  dioxide  which  is 
assumed  to  cover  the  spectral  region  2050  to  2400  em“*  (4. 17  -  4.88  micro¬ 
meters).  The  temperature  range  that  is  considered  varies  from  300°  to 
2100°K.  The  objective  of  the  reported  program  was  to  develop  a  computerized 
program  for  predicting  radiant  energy  emissions  which  could  be  readily  inte¬ 
grated  into  a  flow  field  calculation.  A  description  is  given  of  both  the  radiation 
model  and  the  flow  field  model.  The  described  program  provides  both  the 
spectral  and  spatial  intensity  distributions  of  the  emitted  radiation. 
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Section  I.  INTRODUCTION 


This  report  presents  the  results  of  the  second  phase  of  a  program*  which 
has  been  undertaken  to  predict  analytically  the  spectral  and  spatial  distribution 
of  radiation  emitted  from  an  inhomogeneous,  nonisothermal  system  of  hot  gases. 
Results  are  reported  for  the  4. 3-micrometer  band  of  gaseous  carbon  dioxide  (COz) 
which  is  assumed  to  cover  the  spectral  region  2050  to  2400  cm"1  (4,17  to  4. 88 
micrometers) .  The  temperature  range  that  is  considered  varies  from  300"  to 
2100°K.  It  was  pointed  out  in  the  first  phase  1 1  ]  of  this  study  that  the  program 
was  initiated  because  of  current  interest  in  being  able  to  predict  the  radiant 
energy  emitted  by  a  particular  system  of  hot  gases  in  narrow  bands  of  the 
Infrared  spectrum.  The  objective  of  the  program  was  to  develop  a  computerised 
computational  scheme  for  predicting  radiant  energy  emissions  which  could  be 
readily  integrated  into  a  flow  field  calculation.  Such  a  program  readily  provides 
both  the  spectral  and  spatial  intensity  distribution  of  the  emitted  radiation. 

In  particular,  this  report  outlines  the  current  status  of  the  program  with 
emphasis  on  modeling  vie  radiation  em'  oed  from  the  hot  gaseous  combustion 
products  of  fixed  wing  jet  aircraft  (r  th  mid-infrared  portion  of  the  spectrum. 

It  is  well  established  that  one  of  the  rnaj  •.  '‘ources  of  radUtion  emitted  from  a 
jet  aircraft  is  the  radiation  emitted  from  the  not  jet  plume.  This  is  particularly 
true  when  the  aircraft  is  viewed  from  *  direction  such  that  the  hot  tail  pipe  is 
not  seen.  Since  jet  fuels  are  largely  hydrocarbons,  approximately  15  percent 
hydrogen  and  85  percent  carbon,  which  are  burned  in  an  excess  of  air  or  oxygen, 
then  there  is  considerable  C02  in  the  hot  exhaust  gases.  The  C02  molecule 
possesses  a  very  strong  radiating  band  centered  around  4. 3  micrometers.  The 
radiation  arising  frem  this  vibration-rotation  band  of  C02  therefore  represents 
the  major  source  of  plume  emission  in  the  mid-infrared  spectral  region. 

The  ultimate  objective  of  the  overall  program  is  to  be  able  to  determine 
the  radiant  energy  available  to  a  remote  sensor  from  a  particular  target  or 
class  of  targets  operating  under  a  specified  set  of  arbitrary  conditions.  This 
overall  problem  then  readily  divides  into  three  distinct  phases:  < l )  flow  field 
description  of  the  hot  gaseous  plume,  (2)  radiation  model  describing  the  emis¬ 
sion  characteristics  of  this  gaseous  environment,  and  (3)  atmospheric  modifi¬ 
cation  to  this  radiated  energy.  The  second  and  thtrd  phases  have  been  previously 
reported  |l-4j  on  and  will  not  lie  considered  here  in  any  depth. 


'■•The  results  of  the  initial  phase  of  this  work  are  reported  in  a  previous 
work  (1 1,  which  discusses  the  first  phase  of  the  study  leading  to  the  develop¬ 
ment  of  a  general  spectral  cmlssivity  model  for  carbon  dioxide. 


This  discussion  will  be  largely  conceited  with  that  portion  of  the  program 
which  describes  the  method  of  coupling  the  flow  field  description  of  the  plume 
with  the  radiation  model,  in  particular,  the  flow  field  description,  which  pro¬ 
vides  the  temperature  and  concentration  profiles  for  the  radiating  specie,  simply 
serves  as  an  input  to  the  radiation  model  which  in  turn  is  capable  of  describing 
the  spectral  and  spatial  emission  characteristics  of  the  gaseous  jet  plume.  An 
outline  of  the  approach  that  has  been  taken  is  shown  in  Figure  1. 
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Section  II.  JET  FUELS 


Jet  fuels  consist  primarily  of  hydrocarbons  with  only  small  amounts  of 
other  elements  or  compounds  such  as  nitrogen,  oxygen,  sulphur  and  water. 

These  non-hydrocarbons  are  generally  controlled  very  carefully  to  rather  exacting 
specifications  for  tho  various  commercial  types  of  jet  propulsion  fuels.  The 
total  amount  of  these  non-hydrocarbons  generally  does  not  exceed  1  percent  of 
the  fuel  composition.  Thus,  for  any  practical  calculation  involving  jet  fuels,  it 
can  be  assumed  that  the  fuels  consist  only  of  carbon  and  hydrogen.  Both  turbo¬ 
prop  and  turbojet  engines  use  the  same  types  of  fuels.  These  propulsion  fuels 
can  be  divided  into  three  main  types  [5 }: 

1 )  Aviation  kerosenes  which  have  boiling  temperatures  in  the  range  of 
140°  to  280°C.  Such  fuels  with  crystalline  temperatures  not  higher 
than  -60°C  are  considered  excellent  jet  fuels.  They  possess  high 
heats  of  combustion,  low  saturated  vapor  pressures,  and  good 
viscosity  characteristics  which  ensure  normal  functioning  of  the 
jet  engine  under  a  rather  wide  range  of  operating  conditions. 

2)  Wide  range  distillate  fuels  which  include  gasoline,  kerosene,  and 
ligroine  fractions  and  have  boiling  temperatures  in  the  range  of 
60°  to  280 °C.  This  type  of  fuel  is  highly  volatile  and  generally 
possesses  a  high  saturated  vapor  pressure.  This  causes  high 
altitude  operating  difficulties  because  of  vaporization  and  boiling. 

3)  Heavy  petroleum  fractions  which  have  a  low  vapor  pressure.  This 
type  of  fuel  is  used  principally  for  supersonic  flight  speeds;  and  in 
naval  aircraft  and  training  aircraft  because  of  its  very  higli  flash 
point. 

Table  I  lists  the  most  common  jet  fuels  used  in  civilian  and  military 
aircraft. 

i 

The  C:H  ratios  listed  in  Table  I  were  estimated  from  the  empirical 
relationship5 

I 

H%=  26  -  15pls  ,  |  (la) 

and 

C%=  100  -  H%  ,  (lb) 

where  p15  is  the  density  of  the  jet  fuel  at  15 8 C  in  g/cm3.  The  normalized 
compositions  were  then  derived  by  assuming  the  total  molecular  weight  to  be 
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100  g/g-mole.  In  general,  the  C:H  ratio  increases  somewhat  with  the  heaviness 
of  fractional  composition  and  as  the  C:H  ratio  increases,  the  amount  of  air 
required  for  complete  combustion  decreases.  The  amount  of  air  required  for 
the  complete  combustion  of  a  quantity  of  fuel  is  easily  calculated,  and  for  most 
jet  fuels,  the  ratio  varies  from  14  to  15  units  of  air  to  one  unit  of  fuel.  This 
can  be  illustrated  by  computing  the  mass  of  air  required  for  the  complete 
combustion  of  a  unit  mass  of  fuel.  One  may  consider  a  hydrocarbon  of  the  form 
and  assume  that  air  is  composed  only  of  nitrogen  and  oxygen.  This  is 

reasonable,  since  the  atmosphere  is  composed  of  the  following  major  constituents. 


TABLE  II.  ATMOSPHERIC  COMPOSITION 


Specie 

Mole  Percent 

Nitrogen 

78,09 

Oxygen 

20.95 

Argon 

0.93 

Carbon  dioxide 

0.03 

Hydrogen 

0.01 

This  shows  that  nitrogen  and  oxygen  compose  over  99  percent  of  the 
atmosphere  by  volume.  It  can  therefore  be  assumed  that  the  oxidizer  is  of 
the  form  N^O^.  For  the  temperatures  and  pressures  encountered  in  a  typical 

jet  engine,  nitrogen  acts  only  as  an  inert  specie.  Therefore,  the  following 
burning  or  combustion  process  may  be  considered: 

C  H  +  (X  +  Y/4)Oj  -XC02  +  (Y/2)H20  .  (2) 

x  y 


The  mass  of  oxygen  required  for  complete  combustion  of  a  unit  mass  of  fuel  is 
then  given  by 


Mass  oxygen 
Mass  fuel 


(2X  +  Y/2)Mq 

XMc+  ymh 


(3) 


where  Mj  represents  the  molecular  weight  of  the  i  specie.  In  order  to  com¬ 
pute  the  mass  of  air  required,  then  it  is  only  necessary  to  divide  the  above 
expression  by  the  mass  fraction  of  oxygen  occurring  in  the  atmosphere.  This 
is  readily  obtained  by  multiplying  the  mole  fraction  listed  in  Tabic  II  by  the 
ratio  of  the  molecular  weight  of  oxygen  to  the  mean  molecular  weight  of  the  • 
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atmosphere.  This  gives  23.  2  percent  by  weight  for  oxygen.  Thus  the  mass  of 
air  required  for.  complete  combustion  of  a  unit  mass  of  fuel  is  given  by 


Mass  air _ 1_ 

Mass  fuel  "  0.232 


(2X  +  Y/2)Mq 

XM_  +  YM 
C  H 


(4) 


In  general,  for  the  jet  engine  type  of  problem  being  discussed,  an  excess  air 
factor  is  used  to  limit  the  temperature  of  the  gases  entering  the  turbine.  This 
factor  is  generally  of  the  order  3. 8  to  4.  0.  However,  at  this  point,  the 
oxidizer  fair)  to  fuel  ratio  can  be  simply  denoted  as  O/F,  and  this  can  be 
treated  as  a  parameter.  The  composition  of  the  combustion  products  can  now 
be  obtained  by  first  computing  the  coefficients  A  and  B  for  the  atmospheric  com¬ 
position  N  0  .  As  previously  noted,  the  mass  fraction  of  oxygen  is  0.232  and 
A  13 

can  be  computed  from  the  composition  by  the  expression 


fQ=  0.232 


BM^  +  AM.. 
U  -  N 


(5) 


The  term  "B"  is  uniquely  defined  by  assuming  100  as  the  denominator.  ’  A 
similar  calculation  for  nitrogen  gives  the  result 


N  O  =  N  O,  „ 
A  B  5.48  1.45 


The  composition  of  the  product  constituents  can  now  be  obtained  from  the 
reaction 

(0/  F)  NO  +  C  H  —  XC02+  (Y/2)  H,0 
A  B  x  y  *  i 

+  {(B/2)0/F-X-Y/4)  02+  (A/2)  (0/F)N2.  >  (7) 

If  it  is  now  assumed  that  a  unit  mass  of  fuel  is  being  consumed,  then  the  mass  of 
the  exhaust  products,  W^,  can  be  written  from  the  above  equation,  i.e. , 


(8) 

(9) 


r> 


(10) 


w 


B  (O/F)  -  2X  -  Y/2 


O, 


XMC  +  ™h 


M 


O 


A  (0/F)M 


W 


N 


n2  xmc  +  ymh 


(11) 


For  unit  mass  of  fuel,  a  mass  of  O/F  units  of  air  is  assumed  to  be  consumed  so 
that  a  total  of  (1  +  O/F)  units  is  involved  in  the  reaction.  The  mass  fraction 
of  each  product  specie,  <y  ,  is  then  given  by  simply  dividing  the  exhaust  product 

constituent  mass  by  the  total  mass.  Thus 


W.  W. 

i  l 

“i  ~  2W.  "  1  +  O/F 


(12) 


so  that  the  mass  fraction  of  C02,  for  example,  contained  in  the  exhaust  products 
after  combustion  is  given  by 

x(Mc  * 2Mo) 

®C02  (XMC  +  YMn)(l  +  O/F)  ' 

Since  the  mole  fraction  of  the  combustion  products  is  also  of  concern, 
particularly  if  concentrations  are  expressed  in  partial  pressures,  then  it  will 
be  useful  to  consider  the  relationship  between  mass  fraction,  a.,  and  mole 

fraction,  /?..  For  an  ideal  gas  mixture  in  thermal  equilibrium,  the  total  pressure 

can  be  written  as  where  is  the  total  number  of  moles  and  understood 

to  be  the  weight  of  the  gas  divided  by  the  mean  molecular  weight.  If  the  total 
pressure  of  the  gaseous  misture  is  now  taken  to  be  the  sum  of  the  partial 
pressures  present  in  the  mixture,  then  for  each  constituent  the  partial  pressure 
is  Thus 


or 


Vi 


(14) 


r w.  m 

Pi  =  ptIm!~wI 


M1 
PT°'i  M. 


(15) 
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Equating  these  two  equations  yields  the  desired  relationship,  namely 


(16) 


where  M  is  the  average  molecular  weight  of  the  gaseous  mixture  which  is 
easily  shown  to  be 


M 


_  1 

T  =  Zo’./M. 
i  x 


(17) 


where  the  sum  extends  over  the  number  of  gases  present  in  the  mixture.  Thus 
for  any  specie  within  the  mixtui’e,  the  relation  between  the  mass  fraction  and 
mole  fraction  is 


tti  1 
Mi  2k*. /M. 

l  l 


(18) 


Substitution  of  previously  defined  quantities  into  the  above  expression  yields 


3  .  . 4X  ... 

PC02  Y  +  2(0/F)  (A  +  B) 

(19) 

B  =  2Y 

pH20  Y  +  2 (O/F)  (A  +  B)  ’ 

(20) 

2A (O/F) 

\  Y  +  2  (O/F)  (A  +  B) 

(21) 

2B  (O/F)  -  4X  -  Y 
'  02  "  Y  +  2 (O/F)  (A  +  B) 

(22) 

The  l’esults  of  the  preceding  discussion  can  be  illustrated  by  considering 

a  specific  example.  The  common  hydrocarbon  fuel  CH2  may  be  considered  which 

has  a  C:H  ratio  of  5.96  and  consists  of  14.4  percent  hydrogen  and  85.6  percent 

carbon.  Normalizing  the  molecular  weight  to  a  value  of  100  gives  a  chemical 

formula  of  C  H  .  The  mass  of  oxygen  required  for  complete  combustion 
7.  1«)  14* 

of  a  unit  mass  of  fuel  is  given  by  equation  (3), 


Mass  oxygen 
Mass  fuel 


/l4.  26  +  7.15\ 

\  100  / 


16.0 


3.  43 


The  mass  of  air  required  for  the  combustion  process  is  given  by  equation  (4) 
and  obtained  by  dividing  the  above  number  by  0.  232: 


Mass  air  _  3. 43  _ 

Mass  fuel  0.232 

Using  an  excess  air  factor  of  four  then  gives  an  oxidizer  to  fuel  ratio  of  59.  2. 
Since  this  occurs  as  a  parameter  in  determining  the  concentration  of  the  product 
constituents,  an  O/F  ratio  of  60  will  be  used.  The  mass  fraction  of  each  of 
the  product  species  is  then  given  by  equation  (12) .  Substituting  into  this  set 
of  equations  yields  the  following: 

Product  Mass  Percent 

COz  5. 14 

H20  2. 11 

02  17. 20 

N2  75. 51 


Equation  (17)  can  then  be  used  to  compute  the  average  molecular  weight 
of  the  product  gas  mixture.  This  yields  MT  =  28.85.  Equations  (16)  or  (19) 
through  (22)  then  give: 

Product  Mole  Percent 

C02  3. 37 

H20  3. 38 

02  15.51 

N2  77. 74 


Finally,  it  will  be  pointed  out  that  the  C02  concentration  present  in  the 
exhaust  products  as  computed  above  is  approximately  100  times  the  normal 
atmospheric  concentration  of  C02.  The  mole  percent  or  percent  by  volume  of 
the  constituents,  when  expressed  as  a  fraction,  also  gives  the  partial  pressure 
of  the  specie  when  the  total  pressure  of  the  gaseous  mixture  is  one  atmosphere 
[equation  (14)  J.  Figure  2  shows  the  mole  percent  of  C02  contained  in  the 
exhaust  products  as  a  function  of  C:H  with  air  to  fuel  ratio  expressed  as  a 
parameter. 


MOLE  PERCENT  OF  CO 


2.8  I - 1 - 1 - 1 - 1 - 1 - 1 - 1 - L- 

5.6  5.7  ~*>5.8  5.9  6.0  6.1  6.2  6.3  6.4 

C:H  RATIO 


FIGURE  2.  MOLE  PERCENT  OF  C02  CONTAINED  IN  EXHAUST  PRODUCTS 
AS  A  FUNCTION  OF  CARBON  TO  HYDROGEN  FUEL  WEIGHT  RATIO  WITH 
OXIDIZER  TO  FUEL  RATIO  AS  A  PARAMETER 
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Section  III.  FLOW  FIELD 


The  flow  field  or  plume  model  must  be  capable  of  generating  the  tem¬ 
perature  and  specie  concentration  within  the  exhaust  plume  as  a  function  of  the 
initial  values  at  the  tailpipe  exit  plane  and  the  ambient  conditions.  Figure  3 
depicts  a  typical  jet  plume  along  with  the  objective  of  the  plume  or  flow  field 
model.  An  axisymmetric  coordinate  system  is  used  to  define  any  point  within 
the  plume.  The  coordinate  X  measures  the  downstream  position  from  the  exit 
plane  of  the  nozzle  or  tail  pipe  and  the  coordinate  R  measures  the  radial  position 
or  distance  normal  to  the  symmetry  axis  within  the  plume.  The  plume  is  divided 
into  two  regions  called  the  core  and  developed  or  mixing  regions.  The  core  of 
the  plume  is  a  region  of  constant  velocity  and  thermodynamic  properties, 
whereas  the  developed  portion  of  the  plume  is  that  region  where  the  plume 
possesses  property  gradients  in  both  the  axial  and  radial  directions.  This 
model  is  a  finite  difference  flow  field  program  which  takes  the  conditions  at  the 
exit  nozzle;  velocity,  pressure,  temperature,  specie  or  constituent  concentra¬ 
tion,  nozzle  size,  and  ambient  or  free  stream  conditions,  and  then  computes  the 
temperature  and  specie  concentration  throughout  the  plume. 

OBJECTIVE:  TO  DETERMINE  THE  PLUME  SPREADING 
CHARACTERISTICS.  IN  PARTICULAR  THE  AXIAL  AND 
RADIAL  TEMPERATURE  -  COMPOSITION  DECAY 
RELATIONSHIPS  FROM  A  KNOWLEDGE  OF  THE  GEOMETRY 
AND  EXIT  CONDITIONS  OF  THE  NOZZLE. 
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The  primary  problem  is  to  consider  the  mixing  of  two  streams  which  are 
taken  to  be  parallel  but  which  possess  different  velocities,  temperatures,  and 
chemical  compositions.  Only  the  basic  features  of  the  physical  problem  will 
be  mentioned  in  this  limited  discussion.  These  features  can  be  outlined  by 
considering  a  fully  expanded  axi symmetric  jet  exhausting  into  a  uniform  unbounded 
parallel  stream.  This  is  shown  in  Figure  4  where  it  is  assumed  that  the  pressures 
of  the  two  streams  are  the  same.  The  transport  mechanisms  that  take  place 
between  the  two  streams  include  the  transfer  of  momentum,  the  conduction  of 
heat,  and  the  diffusion  of  species  causing  the  two  streams  to  mix.  In  the  region 
downstream  of  the  jet  exit  plane,  three  distinct  regions  exist.  In  regions  I  and 
II,  the  two  streams  retain  their  original  properties.  Region  III  contains  a  mix¬ 
ture  of  the  two  streams  resulting  in  nonuniform  temperature,  velocity,  and 
specie  concentration  profiles.  The  mixing  process  tends  to  reduce  the  difference 
in  the  properties  of  the  two  streams.  The  flow  itself  can  be  either  laminar  or 
tux-bulent,  and  the  analysis  of  either  can  be  accomplished  by  considering  only 
changes  in  the  mathematical  description  of  the  expressions  for  the  transport 
coefficients. 


FIGURE  4.  DISTINCT  PLUME  FEATURES 


For  the  systems  under  consideration,  both  the  pressure  and  temperature 
arc  relatively  low.  This  results  in  the  reaction  time  being  much  longer  than 
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the  flow  time  so  that  mixing  occurs  before  any  reaction  takes  place  |6).  This 
frozen  type  flow  can  be  analytically  described  by  utilizing  appropriate  transport 
coefficients,  an  equation  of  state,  expressions  describing  the  temperature 
dependence  of  the  species  specific  heat  and  enthalpy,  and  the  conservation 
equations  for  continuity,  energy,  momentum  and  species.  A  lack  of  basic 
knowledge  concerning  the  turbulent  transport  properties  forces  one  to  rely 
heavily  on  experimental  data  for  the  eddy  viscosity,  eddy  conductivity  (Prandtl 
number)  and  diffusion  coefficient  (Lewis  number).  The  most  difficult  problem 
was  the  adoption  of  a  satisfactory  eddy  viscosity  model  to  use  in  the  basic 
numerical  program.  A  considerable  amount  of  effort  was  expended  to  achieve 
a  satisfactory  viscosity  model  for  incorporation  into  the  flow  program. 

Typical  results  from  the  program  are  shown  in  Figures  f>  through  9. 

Figure  5  shows  the  axial  values  of  the  temperature  and  velocity  decay  for  a 
particular  jet  engine  compared  with  the  manufacturer's  data.  This  is  for  a 
compressible  axially  symmetric  free  turbulent  jet  exhausting  into  quiescent 
air  at  sea  level  conditions.  Results  for  an  inflight  aircraft  are  shown  in 
Figures  6  through  9.  Figure  6  shows  the  decay  of  the  temperature  and  C02 
concentration  along  the  symmetry  axis.  This  aircraft  has  a  turbojet  engine 
and  is  assumed  to  be  flying  at  sea  level  conditions  with  a  speed  of  280  meters 
per  second;  The  top  curve  is  a  normalized  plot  of  how  the  temperature  varies 
with  downstream  distance  along  the  symmetry  axis.  The  ordinate  is  the  ratio 
of  the  axial  temperature  to  the  temperature  at  the  tailpipe  exit.  At  a  downstream 
distance  Of  35  feet  the  temperature  has  decayed  by  a  factor  of  approximately 
2. 5.  The  dotted  line  shows  the  ambient  temperature  ratio  and  at  35  feet  down¬ 
stream  the  axial  temperature  is  25  to  30  percent  above  ambient  temperature. 

The  lower  curve  in  Figure  6  shows  a  normalized  plot  of  how  the  C02  concentration 
varies  with  downstream  distance  along  the  symmetry  axis.  The  ordinate  is  the 
ratio  of  the  C02  concentration  along  the  symmetry  axis  to  the  concentration  of  the 
tailpipe  exit.  This  represents  a  concentration  which  is  105  times  the  atmos¬ 
pheric  C02  concentration  (0. 33  percent  by  volume) .  At  a  distance  of  35  feet 
downstream  the  concentration  has  decayed  to  within  10  to  15  percent  of  the 
ambient  value.  Figures  7  and  8  show  the  radial  distribution  of  the  temperature 
and  C02  concentration  at  selected  axial  positions.  These  clearly  depict  the 
diminishing  core  and  appearance  of  the  mixing  region  as  a  function  of  downstream 
distance.  The  ordinate  in  both  plots  is  the  ratio  of  the  plume  property  to  the 
free  stream  property.  Figure  9  shows  particular  isotherms  and  the  core  region 
(shaded)  for  the  same  plume. 


*This  basic  program- was  supplied  to  the  author  by  Dr.  A.  Ferri  and 
Dr.  P.  Baronti  of  Advanced  Technology  Laboratories,  Inc.,  400  Jericho 
Turnpike,  Jericho,  New  Yoijk  il753. 
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FIGURE  5.  COMPARISON  OF  PLUME  FLOW  FIELD  ANALYTICAL  MODEL 
WITH  MANUFACTURER'S  DATA  FORA  PARTICULAR  JET  ENGINE 
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FIGURE  (5.  AXIAL  DECAY  OF  TEMPERATURE  AND  C02  CONCENTRATION 
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FIGURE  7.  RATIO  OF  PLUME  TO  AMBIENT  TEMPERATURE  VERSUS 
RADIAL  POSITION  WITH  DOWNSTREAM  POSITION  AS  A  PARAMETER 


FIGURE  8.  RATIO  OF  PLUME  TO  AMBIENT  COz  CONCENTRATION  VERSUS 
RADIAL  POSITION  WITH  DOWNSTREAM  POSITION  AS  A  PARAMETER 
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FIGURE  9.  ISOTHERMS  FOR  AN  INFLIGHT  PLUME 


The  Advanced  Technology  Laboratories  generalized  frozen  mixing 
program  [6]  is  written  in  Fortran  IV  and  can  be  used  on  any  large  scale  com¬ 
puting  system.  The  input  instructions,  a  listing  of  the  program,  a  typical  set 
of  input  data  and  a  sample  of  the  subsequent  output  are  included  in  the  appendix. 
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Section  IV.  RADIATION  CALCULATIONS 


Alter  the  flow  field  calculations  have  been  made  for  a  particular  aircraft, 
the  radiation  calculations  may  be  started.  The  mathematical  details  and  descrip¬ 
tion  of  the  radiation  model  have  previously  been  reported  on  1 1  ].  Figure  10 
shows  the  geometry  which  is  used  for  making  the  l’adiation  calculations.  The 
plume  is  depicted  simply  as  a  truncated  ice  cream  cone.  The  outer  edge  of 
the  plume  is  defined  by  the  isotherm  T/T  =  1.05.  This  was  found  necessary, 
particularly  for  a  jet  exhausting  into  quiescent  air,  because  of  the  slow  exponen¬ 
tial  of  the  plume  width  in  the  radial  direction. 

A  direction  is  then  chosen  in  which  it  is  desired  to  make  the  l’adiation 
calculations.  This  direction  is  indicated  by  O'  and  called  the  aspect  angle.  A 
plane  is  then  set  up  to  cut  the  plume  in  this  direction.  The  aspect  angle  is  the 
angle  between  the  symmetry  axis  and  the  direction  that  is  chosen  to  make  the 
radiation  calculations.  Zero  degrees  is  the  direction  indicated  by  looking 
downstream  along  symmetry  axis,  90  degrees  is  looking  normal  to  the 
symmetry  axis,  and  180  degrees  is  looking  upstream  along  the  symmetry  axis. 
Within  this  plane,  various  parallel  lines  of  sight  are  constructed  for  actually 
making  the  radiation  calculations.  The  coordinate  Z0  measures  the  distance 
of  the  line  of  sight  above  a  line  in  the  plane  which  intersects  the  symmetry 
axis.  The  distance  measured  along  a  particular  line  of  sight  is  indicated  as 
( .  The  <  =  0  43o«it  is  where  the  line  of  sight  first  intersects  the  plume  surface 
The  coordinate  Xy  measures  the  downstream  position  where  the  calculation  is 
being  made  and  R,,  is  the  radius  of  the  plume  at  this  point. 

.  " 

<r 

'Next  the  coordinates  (X,  It)  are  computed  as  a  function  off  ,  Z0,  a,  Xc 
and  Ry  along  the  line  of  sight. 

Z|  +  [f  sin  cv  -  n/Rq2  -  Zy 
X  =  X„  •+  f  COS  O' 

The  increment  A(  which  divides  the  plume  thickness  into  small  segments  is  auto¬ 
matically  calculated  by  the  computer.  The  flow  field  data  have  been  stored  in 
the  computer  previous  to  calculating  the  coordinates  (X,  R)  along  the  line  of 
sight.  A  two-dimensional  interpolation  scheme  is  then  used  to  calculate  the 
plume  properties  along  the  line  of  sight.  There  is  now  a  dimensional  line 
segment,  the  length  of  which  represents  the  plume  thickness,  with  the  tempera¬ 
ture  and  concentration  of  the  radiating  specie  specified  along  the  line.  The 
molecular  band  parameters  are  now  computed  at  each  specified  point  along  the 
line  segment.  The  emission  calculations  can  now  be  made. 


OUTER  EDGE  PLUHE  DEFINED 
BY  T/T  *  1 .05 


FIGURE  10.  GEOMETRY  FOR  RADIATION  CALCULATIONS 


It  is  not  an  easy  task  to  calculate  accurately  the  radiation  emitted  from 
a  homogeneous  gas.  The  problem  is  obviously  more  complicated  when  dealing 
with  an  inhomogeneous  gas.  A  modified  form  of  the  statistical  band  model 
has  been  chosen  for  making  the  emissivity  calculations.  For  the  inhomogeneous 
calculations  this  was  further  modified  by  using  a  method  similar  to  the  Curtis- 
Godson  method.  This  method  takes  inhomogeneous  gas  which  has  a  certain 
transmissivity  and  replaces  it  with  a  homogeneous  gas  which  has  the  same 
transmissivity. 

In  this  discussion  only  the  results  of  the  method  will  be  outlined.  The 
radiating  band  structure  exhibited  by  the  hot  plume  exhaust  gases  is  actually 
composed  of  many  closely  spaced  or  overlapping  spectral  lines.  In  general, 
to  calculated  the  emissivity  of  transmissivity  of  a  gas,  three  parameters  must 
be  known  for  the  absorbing  or  emitting  molecule.  The  molecular  band  parameters 
are  dependent  on  spectral  location  as  well  as  the  state  of  the  gas.  These  three 
parameters  are  denoted  as  S,  d,  and  y  where  S  is  the  average  intensity  of  a 
spectral  line,  d  is  the  spacing  between  lines,  and  y  is  the  line  half-width.  For 
the  inhomogeneous  calculations,  these  parameters  are  combined  in  a  particular 
manner  and  an  equivalent  set  of  parameters  is  computed.  These  equivalent 
parameters  are  denoted  by  fx  and  v  and  shown  below. 


■’=/Y<'||^1'[Pe'T<,,]dYK 

Y(i)  s[A,T«)]dY(() 


These  parameters  are  computed  in  terms  of  the  basic  band  parameters  and  a 
quantity  denoted  as  Y  which  is  referred  to  as  the  reduced  optical  path  or 
equivalent  gas  thickness.  The  reduced  optical  path  is  defined  as 

Y  (f  )  =  /  — df  , 

0  e(i) 


where  (  represents  the  actual  thickness  of  gas  being  considered,  P 
the  partial  pressure  of  the  absorbing  gas,  and  Pe  is  the  equivalent 
The  equivalent  pressure  is  defined  by  the  relationship 
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where  P  is  the  total  pressure  of  the  gaseous  system,  and  b  is  a  spectral 
broadening  coefficient  which  depends  on  the  overall  gaseous  composition. 


The  transmission  of  the  plume  along  the  chosen  line  of  sight  can  now 
be  computed  as  a  function  of  plume  thickness  and  spectral  location  A: 


r(A,f )  =  exp 


This  calculation  is  necessary  since  the  radiation  emitted  from  the  interior  of 
the  plume  must  be  transmitted  through  a  portion  of  the  plume.  At  this  point 
the  computer  has  available  at  each  spectral  point  of  interest  the  transmission 
along  the  line  of  sight  t .  The  plume  thickness  is  now  divided  into  a  number  of 
small  segments  of  thickness  Ai  ^  (Figure  11).  Average  plume  properties  are 

assigned  to  each  segment  so  that  each  segment  can  be  considered  as  a  homo¬ 
geneous  gas.  The  spectral  emissivity  is  calculated  for  each  homogeneous 
segment.  This  is  computed  from  the  equation 


where  all  the  quantities  on  the  right  hand  side  of  the  equation  have  previously 
been  defined.  The  radiation  emitted  from  a  particular  segment  is  then  obtained 
by  multiplying  the  emissivity  by  Planck's  black  body  spectral  emission  function. 
This  is  then  multiplied  by  the  transmission  factor  to  obtain  the  radiation  emitted 
from  the  surface  of  the  plume  by  the  segment  Af ..  The  result  is  indicated  by 
Ij(A,  0)  and  given  by  the  relationship  1 

e.CjA~5T. 

i  1  i  j 

7r  £exp(C2/AT)-  lj 

This  is  repeated  for  all  segments  along  the  line  of  sight,  and  the  results  are 
then  summed  to  obtain  the  spectral  radiation  from  this  one  point  on  the  plume 
at  the  spectral  point  A.  This  procedure  is  repeated  for  each  spectral  point 
to  get  the  spectral  distribution  of  radiation,  again  from  this  particular  point 
on  the  plume,  in  the  direction  of  a.  The  whole  set  of  computations  is  now 
repeated  for  several  points  on  the  plume  surface  in  order  to  obtain  the  spatial 
distribution  or  the  total  radiation  emitted  from  the  plume  in  the  direction  of 
the  aspect  angle  a: 

KA.cv.Xp,  Zg)  -  ]Tl  [\V/  (cm2-ster-/um)J 
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FIGURE  11.  EMISSION  MODEL 


Figure  12  shows  a  comparison  of  the  results  of  the  homogeneous  and 
non-homogeneous  calculations  for  a  particular  homogeneous  gas.  This  repre¬ 
sents  the  spectral  emissivity  of  a  slice  of  gas  taken  at  the  exit  plane  of  a  typical 
exhaust  nozzle. 

Results  from  the  program  arc  shown  in  Figure  13.  This  is  for  an  inflight 
aircraft  and  for  the  same  aircraft  operating  in  a  tied-down  or  static  position. 

Both  curves  are  for  seal  level  conditions  and  each  is  for  a  90-degree  aspect 
angle.  Again  this  is  for  a  typical  jet  aircraft.  The  emitted  radiation  has  been 
intogradcd  over  the  surface  of  the  plume  so  each  curve  gives  the  spectral 
distribution  of  the  plume  emission. 

it  is  interesting  to  note  that  the  radiation  emitted  from  the  aircraft  in  a 
static  position  is  larger  than  that  emitted  from  the  inflight  aircraft.  The  280- 
meter  per  second  curve  has  a  peak  value  of  1770  watts  per  steradian-micrometer 
whereas  the  static  curve  reaches  a  peak  of  2800  watts  per  steradian-micrometer. 
This  is  apparently  due  to  the  free  stream  air  which  cools  and  reduces  the  size 
of  the  plume.  1'his  is  also  indicated  by  considering  the  area  under  each  curve 
vliH1  represents  the  total  emitted  radiation.  These  values  are  indicated  by 
Figure  13.  where  it  can  be  seen  that  the  static  value  is  approximately  30  percent 
less  than  the  inflight  value. 
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FIGURE  12.  COMPARISON  OF  HOMOGENEOUS  AND  NON-HOMOGENEOUS 
CALCULATIONS 


The  absorption  dips  occurring  in  each  curve  are  due  to  the  cooler 
outer  extremities  of  the  plume  absorbing  the  radiation  from  the  hotter  interior. 
This  effect  is  more  pronounced  for  the  static  curve  since  the  plume  is  larger 
in  size.  The  small  peaks  seen  beyond  4.  8  micrometers  are  due  to  a  much  weaker 
C02  band  centered  in  this  region. 

Figure  14  shows  data  for  the  same  inflight  aircraft  that  was  discussed 
in  Figure  13.  Again  this  is  for  a  90-degree  aspect.  The  results  are  presented 
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FIGURE  13.  TYPICAL  TURBOJET  PLUME  SURFACE 
EMISSION  AT  90-DEGREE  ASPECT 


in  a  different  manner,  in  that  the  spatial  distribution  was  not  integrated  over  the 
whole  surface  of  the  plume  but  only  over  the  radial  coordinate  and  then  integrated 
over  the  spectral  region.  The  result  gives  the  axial  decay  of  the  radiation  or 
the  radiation  per  unit  length  of  plume.  This  permits  the  radiation  centroid  to  be 
calculated  which  is  also  shown  in  Figure  14.  The  plume  geometrical  centroid 
is  about  7  feet  aft  of  the  tailpipe  for  this  particular  aircraft. 
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FIGURE  14.  AXIAL  DECAY  OF  PLUME  RADIATION  FOR  A  TYPICAL 
TURBOJET 


One  of  the  major  objectives  of  this  effort  is  to  be  able  to  predict  the 
radiant  energy  that  is  emitted  from  a  jet  aircraft  and  is  available  to  a  remote 
sensor.  This  implies  that  the  signature  must  be.propagated  through  the 
atmosphere  which,  in  reality,  aets  as  an  attenuation  filter.  Therefore  an 
.  atmospheric  transmission  model  (3]  has  also  been  developed  in  order  to 
determine  the  atmospheric  modification  to  the  emitted  energy.  It  is  a  two- 
component  model  that  was  derived  from  the  well-known  experimental  data  of 
Taylor  and  Yates.  Molecular  absorption  of  HjO  was  considered  along  with 
other  constituents  of  constant  concentration  such  as  C02,  N20,  CR,  and  CO. 

The  constituents  of  constant  concentration  were  combined  into  a  single  species, 
X02,  The  amount  of  X02  was  chosen  to  be  32  atm  cm/km,  which  is  the  average 
concentration  of  C02.  Figure  15  shows  a  comparison  of  the  transmission  model 
with  a  set  of  independent  measured  transmission  data.  It  can  be  seen  that 
the  computer  model  matches  the  measured  data  quite  nicely.  The  broad 
absorption  dip  between  4. 2  and  4. 45  micrometers  is  due  to  C02  absorption.  The 
dip  at  4. 5  micrometers  isduetoNjO,  and  the  fine  structure  is  due  primarily  to 
HjO. 
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FIGURE  15.  COMPARISON  OF  ATMOSPHERIC  TRANSMISSION 
MODEL  WITH  MEASURED  DATA 


Figure  16  shows  a  typical  plume  signature  put  through  1,  7,  and  15 
kilometers  of  atmosphere  by  using  the  previously  described  transmission  model. 
This  shows  the  importance  of  the  wings  or  skirts  of  the  emission  band  since  the 
center  of  the  band  gets  completely  absorbed  in  relatively  short  path  lengths. 

A  comparison  of  the  results  of  the  present  model  with  measured  data 
is  shown  in  Figure  17.  This  is  for  a  particular  jet  aircraft  operating  in  a 
static  position  at  100-percent  power.  The  aspect  angle  is  10  degrees  and  the 
measured  data  were  recorded  at  1  mile  from  the  aircraft.  The,  theoretical 
results  were  put  through  the  above  described  atmospheric  transmission  model 
to  obtain  the  signature  shown  in  Figure  17. 
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FIGURE  16.  REMOTE  PLUME  SPECTRAL  SIGNATURE 
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FIGURE  17.  COMPARISON  OF  MODEL  WITH  EXPERIMENTAL  DATA 
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Section  V.  RADIATION  COMPUTER  PROGRAM 


A  copy  of  the  computer  program  is  shown  in  Table  III.  A  typical  set  of 
input  data  is  shown  in  Table  IV.  The  first  72  data  cards  are  standard  input 
for  each  run.  These  cards  alternately  list  the  molecular  C02  band  parameters 
'A 

S/d  and  S  2/d  as  a  function  of  wave  number  and  temperature.  The  first  data 
record  on  each  card  gives  the  wave  number  and  the  remaining  seven  records 
list  the  band  parameter  at  the  given  wave  number  as  a  function  of  temperature 
in  increments  of  3006K  (for  second  data  record,  T  =  300  °K;  for  third,  T  =  600°K; 
for  eighth,  T  =  2100eK). 

The  next  two  data  cards  are  also  standard  input  for  each  run  and  simply 
contain  headings  (RADIATION)  for  the  output  data  headings. 

The  next  two  cards  are  used  to  document  the  date  of  the  run  and  the 
description  or  title  of  the  run.  The  first  card  contains  the  date  (first  12 
columns)  and  the  second  the  title  or  name  of  the  run  (first  72  columns).  The. 
next  card  must  contain  two  pieces  of  information  according  to  the  format  (215). 
The  first  is  the  number  of  X-stations  (NXS)  that  are  desired  in  the  calculations 
and  the  second  is  the  number  of  radial  or  Z  points  (NZPEX)  desired  at  each 
X  station.  (These  numbers  must  be  less  than  or  equal  to  60  and  8,  respectively). 
The  computer  checks  to  insure  that  these  limits  are  maintained.  The  coordinate 
X  measures  downstream  position  in  centimeters,  and  Z  is  the  coordinate  normal 
to  X,  also  measured  Tn  centimeters.  The  actual  magnitude  of  these  numbers 
and  their  spacing  is  determined  by  the  flow  field  computation  Again  NXS  and 
NZPEX  count  only  the  number  of  the  points  at  which  calculations  are  made.  The 
set  of  flow  field  data  is  next  read  into  the  computer.  This  is  denoted  between 
statements  172  and  110  of  the  program  listing.  An  axisymmetric  coordinate 
system  is  used  so  that  any  plume  property  is  defined  by  the  coordinates  X  and 
R.  Data  from  the  flow  field  calculations  are  written  on  tape  at  specific  X  or 
downstream  locations.  The  number  of  radial  points  at  each  X  position  is 
denoted  (NRP)  and  each  radial  point  defined  by  R(I).  At  each  downstream 
position  where  a  radiation  calculation  is  desired,. The  temperature  T(I,  J)  (I 
denotes  the  number  of  the  downs  tream  -posi  tion,  J  denotes  the  number  of  the  radial 
position)  in  degrces-Kclvin,  the  radial  position  R(I,  J)  in  centimeters,  and  the 
partial  pressure  of  C02  (PC02)  in  atmosphere  is  written  on  tape.  The  num¬ 
ber  of  radial  points  at  any  given  X  position  must  be  less  than  or  equal  to  25.  It  is 
assumed  that  the  total  pressure  anywhere  in  the  plume  is  1  atmosphere.  These 
quantities  are  not  written  or  read  according  to  aformatbut  rather  in  binary  notation. 

Finally  a  card  is  read  which  contains  two  quantities,  ALPHA  and  ISTOP 
according  to  the  format  (F10.  5,  15);  ALPHA  is  the  aspect  angle  in  degrees  and 
ISTOP  is  either  1, 2  or  3.  A  (1)  will  terminate  the  run  after  the  calculations 
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TABLE  III.  RADIATION  COMPUTER  PROGRAM 


-  CFN  SCLRCE  STATEMENT  -  IF  MS  F 


TRACY  JACKSCN  7/24/68 _ _  _ _ 

CO  2  A. 3  M1CRCN  BAND  PLUME  SPECTRAL  RADI  AT  I CN  PROGRAM 

Cl  ME  NS  ION  .X.(tC).R|60,2S),T(tC.2;).,PCC2l6O.25).SOCI36,7).SOU2l  36.7) 

UI  PENSION  WAVNC I  36 )  .VtAVLHI  36 1 •  TT  1 7 1  <  TSPRAD 1 36 ) 

READ  HAND  PARAMETERS  FCK  C02  EM  I  S  £  I V  I TY  CALCUIATICNS  — 


REMIND  10 

CO  1M  U1.36 _  _ _ _ _ 

REAP  15  •  1000  )  MAVNOII  ),  I  SCO! I, J)  ,J«1.7) 


ICC  mAVLHU)*10CCC./XAVN011) 

JOJiO  FORM  AT  I F7. 1 .7ES. 2) _ 

1X1  1C1  j«1.7 

_TJ«J _ _ 

TT  C J  I*  300  •*  T  J 

_ no  l ci  i»l.3t _ 

1C1  SHC2II  .J)«S002II.JI*SCO?ll  .J) 

DIM;; NLIliN  DATEI2).TI  TLE!12  I HEACNC 1 1 H I 
Kb  AT)"(  S.  10 5 1 )  I hEADNGI  Jt  ,J»1,18) 


1  form  at  1 12 At /t At) 

RE  AD  IS. 10S0  )  CATEI l ) .DATE  1 2 ) . IT  I  TIE  I J ) . J«1 .12) 

AC 

»R  1 T  E  ( 6.1  IOC  I  OATE(1),OATE(2».  1  T  ITLE  (I ) ,  I«  1. 12 ) 

11  CO  FORMAT!  1H1.ACX.52HCU2  A. 3  M1CKUN  EANC  PLUME  SPECTRAL  RAC1AT10N  PRO 
‘lGPAM7lH0.42X.2A6/lFC.3CX.l2At) 

MR  IT  El  4 . 1 1C  1 1 _ _ _ 

TlCl  FORMAT UHC.A1HALL  FLUME  C0CR01NATES  a‘rE  MEASURED  IN  CMr) 

_ FRjT.LLt.tU.qil _ 

11C?' FORMAT ( IHO. 17FMAVENUMBER-I 1/CM) ) 

_»R  IT  El  6  ♦  1 10  3 ) _ 

TlC3  FORMAT!  IHO , 2CF.MA VE LENGTF-I  MICRONS  1 1 

_ HR1TEU.110A)  _ _ _ 

11CA  FORM  AT  1 1HC.  AIFRAOI  ATICJN-IMAITS/SC  .CM^-MICRCN-STEPAOI  AN ) ) 

NXS«NU“HFR  X  STATIONS 

■_  N/PFX^NUHBCR  2  PC1MS  AT  EACH  X  STATIGN _ 

RE  ADIS.  1001)  NXj»N2PEX _ 

loci  FORMAT (215) 

_ go  tc  it; _ 

IFIN2PEX.GT.E)  GC  TC  171 

_  Oil  TO  172 _ _ 

1*70  MR  IT  El  6.1 120 

1 1  2C  FU Rf  ATI1H1.A0H  NUMBER  OF  X-STATI C NS_MUS_T  '  BE  60  CR  LESS) _ 

STOP 

171  MRITEIt.1121) _ 

1121  FORMAT  1 1H1,  26FNUMBER'  OF  2-PCINTS  MLST  BE  8  CR  LESS) 

_.sjrp _ _ _ 

Z  _  READ  INPUT  DATA  TAPC  FROM  PLLME  FLCM  FIELD  PROGRAM _ 

Z~"  'MpiNUM.BET  CF  RADIAL  POINTS  . 


172  CONTINUE 

RcADtlC)  I T  ITLEI 1 1  »I»1 . 12) _ 

RcAD 110)  NCRUN 

no  ic?  i» i.nxs _ _ 

RF  AO I1C)  XI  I )  ,NkP 

READ  1 10)1 T (  I  ♦ J ) ,RI 1 . J) . PC02 1 1 . J ) ,  J* 1 .NRP) 
IF  INRP.LT.2EI  GO  TC  110 

_ GU  TC  102 _ 

110  NRP1*NRP+1 

Oil  1  11  KMNRF1.25 


A 


0000  0)0  j  ooooo 


TABLE  III.  RADIATION  COMPUTER  PROGRAM  (Continued) 


TtST  -  f.fN  SCUWCE  STATEMENT  -  CFMS )  - 


m  I.KK  »  =R  (  I  ,NfP) 

Tl  I.KK  1  =  1 1  I  ,NPP)  _  _  _ 

111  PC. 21  1  ,Kk  )  =  PCC2I  t.  K-WPJ 

l_C?  (.  iNT  1MUC _ 

.  IPAGE=G 

AL  PliA*  ASPEC I  ANOLc  IN  UEGKEESCO  CtG*NCSE-ON) 

IS  llfPs  I.aCLL  ;  EKN»I  NATfc-  PROGRAM  AFTER  SINGLE  ALPHA  CALCULATION 

IS TllP®  2.KET  URNS  til  THIS  POINT  FCF  ANCTHER  ALPHA 

ISlilP-1. RETURNS  rn  STPT  173  TO  PEAL  NEW  SET  PLOW  FIELU  DATA 


UC 


ni. 

ill 


ixr=  i 

Ri>  *f  1 1  Xl)«  2S 1  _  _  _ _ _ _ , 

xij  *x  ( t  xTi  >  t  •  c  1 

If (ALPHA. PU.tEC. Cl  XG-XlNXSl _ _ - _ 

ZPf»*N/PFX 

Ilf  U  MC/2PE  X _ _ _ 

ilil  6C2  »=l.NZPtX 

2M  =M  _ _ _ _ _ _ _ 

/)»nEL/*('fn-i.l 

Ip  I  ALP  HA.  C  w‘ .0 .0 1  RC*2.»KCI  »?5>  ♦  1/N-l .  )  »( K(  NXS.25  »-2.«ft(  1.25)  > /ZPEX 
lHALPHA.EO.lEO.CI  KC=(2P-1.I*(H(NXS.25))/  ZPEX 


1C3  'EAIHS.IoICIALPHA.  I  S TOP  _ _ 

lulG  fORy.ATIF10.Si  151 

AL  F»  Al  PHA/  57. 29*>  7HC  _ 

CC*C0STALP|  "  . *  -  . ' . . 

SC»S1N(A1  PI _ 

Fi.i.i«is92e 

CM  It)  L«l’.  3t  _ _ _ 

i  to  rsPKA.)(u'*o. 

NXSM1-NXS-1 _ _ _  _ _ 

if  (ALPHA. EO. 6. C.CR. ALPHA. cc.-lao. Cl  NXSP'l*! 
C  I  1  06  |c  1.NXSP1 _ . _ 


........ ....*.»**♦».»«..»*.♦****..»********•»*•********< 

t\ LOLL  ATP  CCCH01 NA TfS  X.R  ALCNG  l INE-CF-SICHT  AS  FUNCt ION  OP  ITS 
LENGTH  -  -  -  -  -  •  • 


oiypNsiLN  cr(  ioyi.cxitcci.CLCfoci 

CELL  »IC. _ 

If  (ALPHA.  10.  C.C.C’t.  ALPHA.  EC.lbC.  Cl  DEL L*X I NXSI/08. 

2CC  XL*0. _ _ _ _  _ _ _ 

M  2 G 1  Jal.iOC "* 

_  _ _ _ 

APf  =  XL*SC-SGPf  (KC«RC-Z0»Z.1I 

_  CP  CN  1  =  SCHTC  ZC« 7C»APP»AFR| _ 

C* IN  1-XCfXL  *CC  " 

•  Cl.  IN )*XL  _  _  _ _ 

CMCXCM.LT.C.)  GU  TO  2 CA 

If  1C  X(  M.CT  .XCNXS)  I  GU  TC  2CA__  _  _ _ 

If  (ALPHA.'fc'J.O.O.GR. ALPHA. EC. 180.CI  GC  Tl)  203 

CO  2  Cl  K»1,NXS _ ....  _ _ _ 

Kl  *K 

I  flCXCM.Lf  .XIMI  GC  TC  202 _ 

201  CuNTINUE 

202  |F(f,PIM.GT.R(Kl.?SI>  GC  TC  2GA  .  ...  .  _ 

20  XL  -XL*OfcLL 

C  _  _ 

C  DILI  T i,()  SMALL, INCPfASt  AS  fClLCVS 

C  _  _ 

CeLI  s')tLl  f'JLLL/2. 

GU  TC  2uO  _ 

2CA  If  CM. LT. SOI  GC  TC  2CS 

G'  IC  200  _ 


TaF 

IAS 
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TABLE  III.  RADIATION  COMPUTER  PROGRAM  (Continued) 


-  bbl.  SOURCE  STATEMENT  -  IF  M  S ) 


_ J.OLlA&Biji£Kfet*5 C  A5_LaJL£*L 5 . 

2C5  CE  LL  *C L  IN  I  /  7 E  • 


PLLilfc  THICKNESS  NEGLlGlHLE  At  IHS  PCIM.MCVE  UN  TQ  NEXT  POINT 
AFTER  SETTING  ALL  SPECTRAL  RAG  I  AT  ILN  VALUES  EuUAL  10  2ERO 


IEIuELL.Lb.C.CCS)  GO  IC  41C 
GO  TO  2  JO 


C*<  **<■*>.  ♦*  *«  *»*<<»****,***»»»»**#*i.*,.*»*»*********»*****»****»t**  ******** 

C  i.ALCULATE  PRESSURE  AM  TfcMPEKA T LE t  ALONG  l INE-OF-S IGHT ' AT  EACH  X,R 

cYrV asY>Tn~c tTTcoi Tc p i’i cfi  ~ 

2Q9  00  3C0  J«1.N 
GO  JCi  i*l.  NXS 
NN  =1 

IF  K. XI J  i.LE  .X  III  t  or  TC  302 

3C1  CUNT  INGE _ _ _  _ _ 

TC?  CO  3C3  L»l,25 

VM«L _  _  _ _  _ _ 

IF  lu  Rf  JT )  •  LE .HINN »L  j”!  GO  TO  304  . 

303  CONTINUE  _  _ _ _ 

3C4  CKl«Cft(Jl'-K(NN-l,'tP-U 

0k7*R(  NN-l.MHI-RINN-l.PK-t  I  _  _ _ 

nR3=CRIJ»-KlNN,PM-ll 

CE  4*fU_NN«PN  I-RlNNtPP-l  )  _ _ 

cxiicTi  j'i-xTnV-iY" 

CX  2*  XI NM-X I NN-1 1  __J _ 

1xT*'fl  NN-1,  PP-'l>*I  TtNN-1  ,.MP)-i  <  NN-l  .PM-l  1  )*0R1/UR2 

TX2*  TI HN.Prt-1 1*1 T( NN.HPI-T INN. HP- 1 ) J  A0R3/0R4 _ 

CTIJ  )«TXUI  TX2-TX1  »*CX1/T)X2' 

PX  1*_PC02I  I ♦  IPCC2I  NN-1 .  PP  I-PCC2  (NN-i»MM-l) )  *UR  1/DR2 _ 

PX 2*PC()2I  <vN«MP-f  IK PCC2 1 NN  «FP) -FCC2 1 NN.HK-l I ) »0K3/0R4 

CP  (J  )-«PXl*l  PX2-PX1  )»UX  1/0X2  _ _ _ 

iriCTlJ).LT.3CC.I  CT  IJ  j *30C* 

_ 3  CO  lEICPIJI.LT.  C.C00331  C  Pi  J )  *G.OJC33 _ 

“c  " 

0  ***********************  ************************************ ************ 

CALCULATE  AT  f ACH~ SPfcC TRAL  "FcTnT CCN'sfCERtir'f hTIanD  PARAMETER'S 

C  AT  EACH  X.R  ALUNG  II NE-CF- SIGHT  _ 

C 

CIPliNSIUN  jCS0COc.ig0)_.CS002l36.Ki-l  _  _ _ 

00  4C1  K=  1 .  N 

_ IRI»CTIK)/3CC. _ 

1K1RT.LT.U  1RT  =  1 

CT  IK  C  TIKJ -  T  T  I  JR  T IJ/JOG.  _ _ 

on  4  Cl  L»1.36 

CSCu  (L  «KI_;SCO_IL«_tA  Tl+O  T 1* I  SCO!  L  » iPTt  1  l-SCG  I L  » IRT  IJ _ 

4  0 1  CS002I  L.K  I  =  S0C2 1 L,  f«  T)  ♦ilTl *  I  SCU)2 1 L ,  l  RT*  1 1 -SCC2 1 L ,  IKT )  1 

„C _ 

C  **"  ***  *****»****#**»*************»***«<*»**#**********  ******  •****•*«•*** 

C  CALCULATE^  TRANSMISSION  AJJEACH  .SFECTRAL  POINT  FRCM _EDGE  UF  PtOMb _ 

C  TO  X.K  ALONG  L  INE-CF-S IGHT 

C  _ _ _  _  _  _ _ 

01  PENSION  Y 1 1CC1 « 7  RNSI  36,  IOC  I 


CO  5C1  K*  2,  N 

SCI  Y!  K1  *Y(K-l)*Q.;«tCHKl-CHK-lM  « 1CFI Kl »CP( K- 1 )  ) 

06TC2  K*lV36  . . . 

TftNSIK.UK.C _ _  _ 

p*o. 

C*C. 


o  o  o  o;o  j  '  o  u  o  o  vj  o  u 


TABLE  III.  RADIATION  COMPUTER  PROGRAM  (Continued) 


tE^T  -  ir.\|  Sl.LKCH  STATEHEWT 

n,i  scj  j-?.n  ” 

C¥  =.).5*(YIJ  1-YlJ-l  )l 
PV  I=C.  5*tf  IIJ  l«CII  J-l)  ) 

GAP*  C,L  /5*-$L.R  1 1  3C  )  ,/PVT 1 
r»=p»CY»tcsm'(K.j-i  )*csi:cik,ji  i 
_  C  =  OGY»(.AM«  ((■  SOi) 2 C  K«  J- 1  I *C  SCG2  I  K  ,  J)  )  _  _ 

An  i;*  -2 .  m’>j/  p  t  •  i-  l.  *scr  t  ( i .  ♦  p»p7  c  n 

502  KSSIK,JI=t  xp(ARG) 


!*»•««•,*  «*  <■«»«  •  «J*  »**  *«J*»  ♦*»•**»****  •**»♦»**  ***»*»**>■*«**»***•»*****••* _ 

LA  LC  (it  A  TE~  TOTAL  SPTcfwAL  l P ISS I CN/LNl T-APF A  EROW  THE  POINT 
XG.m\7G  IN  THE  i)I»LCTICA  CF  ALPHA 

301 

0 1  ME  NSlf>  bPRAOl  3h,l  1)  ,/RI  II.KSUf  SI9I  .OFLTALO) _ 

IN'>X  =  0 

Gi!  TG  fc20  _  _  _  _ 

A  10  INGE  X-  100 

6?C  /-X  (>|  1  =  /t; _ _ _ 

1KALPHA. EO. C.C.GR. ALPHA. EC. 180.0  Zfi(P)*KC 

NSTfP$l'M*N  _  _  __  _ 

GELTALIwI-OlLI  '  ’  '  "  “  ‘  ‘ 

PX  =M_  ___  _  _ _ _ 

oo'Vw  k»w  it  . . 

8LC  5PaAClK.Ml»C. _ __________________________ _ 

IE  I  INDEX..", I  .fed  GO  TC  fcC2 

l*P  i*A-  1 _  _  _ _ _ _ 

GO  GCl  K=1,NP1 

PVI»(CTCk)»CT(K»1I  )  =0.  5  _ _ _ _ _ _ 

XXMATtie./PVT 

CD  AC  1  L»1.3t _ 

L-'SV  TY-ITRNSIL.iU-TPNSIL.K+IU/TPNSIL.K) 

Pv-L*  HA  VLM I  L  1  _ _ _ 

PI  KF  All*  1 1  ROE.  **  I  -  5  i  1/  I  EXP  (  XX/FRl  1-1. 1  33A  335 

,  «P.'lAC(L,M)*JiLKKAO»E.lSVTm_AAStL.J‘LtSPRAO<L.K) _ _ 

6ci  cuaunlf 

t>C2  COAT  INGE  _  _ 


PEKFQRH  SPATIAL  INTEGRATION  OVER  FLJPE  SURFACE  TC  OBTAIN  THE  TOTAL 

SPFCTRAL  EMISSION  10  ThtJllJitCT  !LA  CF  ALPHA, ALSU  PERFORK_AN  _  _ 

IN  TF  GR  ATI  ON  OVER  THI  SPECTRAL  REGICN  TO  OeTAIN  TEC  TOTAL  RAUIATIUN 
THAT  IS  EMITTED  FkC*  THE  PLLPE  AI  ASPECT  ANGLE  ALPHA  _ _ 

svx  +  i 


2fil“XMI«K<  1  XL"  ,25) 

1)0  7  CO  1  =  1,  3fc  .  _ _  _ _ 

SP5AG(L.RX*1)=C. 

7 CO  SPP.AIX  L,MX»31«C.  _  _ _  _  _  _ _ _ 

1M  ALPHA.  Ew.o.C.CK  .ALPHA.  EC.  180.0  2B(  HXM  )*RCNXS,25) 


0  I  70  P*«!«PXP1 

7C1  SPnACIl  .HX*3)»SPRAUU.PX*3)  +  UBIH  WHIP-11  1MSPRADT  L.H-1  IhSPKAUCL, 
l*>  1 

iriALPHA.EO.C.C.CR.ALPHA.EC.iaj.L)  GC  TU  7C5  _  _ 

IF(ix-).FO.l)  GL  TO  7C3 

_ _ L?Ja-LE _ _ _ _ 

7C2  TSPn  A.'ML)  =  TSP«Ai)<L  I*  (XI  I  l-Xl  I-l  1  M(SL/2.)*ISPRA|)(L.MX»2)*SPRA0(L,M 

1X«  2)  I  __  _  _ 

70  CO  7CA  L=  1 .  3E 

7 CA  SPRAC(L.MX»7|=SPRAC(L.PX»3)  _ 

GO  TC  7l 7 

7C5  CD  7Ct»  L=l,3t _ _ _  _  _ 

GO  7 C8  v  =  2, PXP 1 

7Ct  ISPt  AO  I  L )  s  T  SPhAOIL  )♦  WmIPI  *2rtl  f  I -7E I  *-l )  *1  C(  P-1  I  )*P  I  *(  SPP  ADI  L  .  M  l  +  S 
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TABLE  III.  RADIATION  COMPUTER  PROGRAM  (Concluded) 


SLURCE  STATEMENT  -  IFN15)  - 


— jjm/u _ 


IPKAOtL.M-l)  1/2. 

7C7  on  7C8  L* 1 . 7t 

7CH  3°KA  Dl  L  «MXP  1 1 *SPRAD t L. MXA3 ) 


c««**v *«*«♦*♦*#***♦ »*»** «a»i *«*«*»*•*  «« «»***•**• *»•*•««*»**»*> *««**••*** 

C  '-.RITE  DESIRED  OUTPUT 
C 

IP  Al»E*  IPAliE  ♦  1 

„K  1 7  E (  A. 1200  I1A  TE  II  I.  OATF 1 2 1 »  I  F  AC  t 

12  00  FORMAT  (  IH1.  t2X.2A6.3CX  .AHPAGE. 1 31  _ 

hR  ITT(6tl2Gl|  !TlTlElU)',j'*l  rixT  . . 

«R  ITE!  6,12H5)_NCKUN 
12«5  format' iVho»6ox»ThkVn  NC*.IAI 

12C1  Form  AT  ( 1HC.  3CX.  12A  fcj _  _ 

AO  XX  *NSrEPSril-l 

80T?«ASXX*0F.LTALtl  ) _  _ _ 

HR  IT  F  ( 6,1 20  2 1  ALPHA, X0.R0T2 

12  C2  K>  HI.!  AT  <  1HU ,  1  AH_A  S£ECJ_  ANOL  E«.E1C.A,5X,  2CHDCh'NSTk  FAR  POSIT!  U'lf_.  E10. 

1A,  Ex  ,16HPLUME  TMICKNtSSV.Fl6.Vr*  . 

HRlTE(0,150ajNSTEJS(J»,J*l.HX)  _  _ 

15C0  FORMAT!  irtO,  7h  STEP $» .A X. 91 12 1 

_ »R  1TF(  e.1501  1 10FLTAL!  J  I  ,J»1  ,PX) _ 

13C1  Format!  1M  .  SH  DELTA  L*  ,  7X  ,9E  12  .  A ) 

_  HR  IT  Lift. 120311 2H(J».J«l.PXl _ 

~120  3  FORMAT !  IH  .17H  PAOIAL  PCSI  T ICM*  .E  1 1 .  A ,  8fc  12  .A  » 

RX2«  2»PXP  1 _ _ _  _ 

wR  1  T~E!  t>»  12GA )  IHEAONGI J  I « J*1  .PX2  1  ~  "  ‘  * . 

12 CA  FORMAT! 1HC.17H  HAVMl&g  »A VLKGH ,2 X < lBAbl _ _ _ 

DO  1205  0«l,3t 

12 C5  >R !T£lt.l20tl  HA VNC! J1 .HAVLH 1 J 1 . 1 SFKAD ( J ,K ) . K»1 . PXP U _ 

1206  FOkMATIIH  .F7.C.F0.3.9E12.A) 

01  MENS  ION  RUK1TL16C1 _  _  _ _ 

RUNI  TL  ( I )  *0  . 

_CU  _125C_J*2._lt _  _  _ _ _ _ 

12  50  P'.i'lfTL  ( 1 1  “RLNITL  (I  !♦  (hAVLH  1  J-l  1  -iTaVLHI  J1 1  *0. 5*!  S FkAO!  J-l ,  MXP 1 1  ♦ 

ISP  RAO!  J.MXP11  > _ 

HR1TF16. 12511  RllM  TL!  I  ) 

12  51  FJRR ATI  IHO.  A6HKAU1  AT  KlN/UNt  T  LENGTH  OF  PLUPFtwATTS/STER-CMl*.  til. A 

“‘*11 

«RITE16. 12071  _ _ _ _ 

12 C 7  FORM Af 1 1H0»  2CX.SAHNCTF-THE  LAST  HAcfATIUN  COLUMN  ABOVE  GIVES  THE  S 

_ 1PFCTRAL  HAOlATlllN/U.HT  LENGTH  OF  FLUFF  SURFACEI _ 

1 06  CONTINUE 

IPAGF.IPAGFaI _  _ _  _ 

HR  It  FI  6, 120  Cl  DATE  ( 1 )  .DATE  12) « 1  FACE 

_ «8  HE(t.l2011ITITLE!0)  ,J»ltl21 _ _ 

«R  ITEIt, 12081  AlPHA 

12C  B  FORMAT!  1H0«  1AH  ASPECT  ANGLE*  .F  1  C  .  A  ) _ 

TI1TKA0*0. 

_  11  1  7C9  L»2. 36 _ _ _ 

7C<T  T'ITp  Ai)«TOTRAOFlHAVlHIL-l>-HAVLHIL  ll*rfSPKAC('i-ll*!TSPRACILIIl 

T'l  TR  A,)»T(|TR  AO/2. _ _  _  _ 

hRIT"E(  t.*12uSl 

12  C5  FORM  AT! IMP. 37P  hAVNHUK  HAVLNGH _ SPtCTKAL  RAC1AT ION) _ 

’  HR  fTE! 611210 IkAVNOI 01  .HAVLHlJI  ,1  SF*f"AD!o1  .  J*1 . 36  1 

1210  FORHATIIH  .F7.C.F9.3.AX.E12.A) _  _ _ 

HR  if  f  I  6. 1 2 11 1*  TUTR  AO 

centr»c._ _ _  _ _ 

nil  1252  I  *l.NXSMl 

12  52  ROM  TL  1 1 1  «RLN  1  TL  1 1  1  *X  1  1 1  »SC*SC _ 

OH  125’  I *2 .NXCM1 

12  53  CeNTP«f.FNTR*(X(I  I;  Xl^- 1 1  Mp. 5* «_RL F  m  1 1  =  1 1  *KUM TL!  m _ _ 

CF  NTK*CEN  TP/TCf  it  All* 

HR  IT  E 1 L .  1 25A )  CENTR  _  _  _  _ _ 

125A  Ft  J  u  V  Af  1 IHO*  SHCEHTRC1U*  «E  11  .A  i 

12  11  FORMAT!  1H0.2AHTJTAL  RADIATION  C  H  T  T£0»  .E 10. A ,  11M_  HM  T  S/STER) _ 

GO  ICC  1C7.  103.172)  .  I  STOP  “  ‘  . . . 

1C7  RJ_hI_N0_10 _ 

ST  rp . * 

CMC 
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TABLE  IV.  TYPICAL  SET  INPUT  DATA 


SOATA 

2050.  2.606-03  4.BCC-Q4.  3.406-04  4.10E-04  6.P0F-04  5.706-03  2.606+02. 

205Q.  A.60C-02  6.60E-02  t. 336-01  3.35E-01  7.90E-0J.  2.90E-00  B..90E-00- “ 

2CA0.  H101E-Q3  1.O3E-03  5.72E-04  6.37E-04  1.37E-03  9.06E+O3  3.56E-02 

2060.  5.25E-02  I.50E-02  I.49E-01  ...766-01  l.OOE-OO  '3.60E-00-  1.00E+01 

2070.  1.87E+02  I.8VS-O1  9.57E-04  1.06E-03  2.28E-03  1 .426-02"  4., 76E-01.' 

.2070.  6.00E-02  8.42E-02  1.68E-01  4.10E-01  1.2BE+00  4.40E-00  .1.13E+01 

2080.  li'5E-0P  1.9OE-03  1.3OE-03  1.66E-03  3.el5;032.13E-02  6.286-02' 
2080.  6.67E-02  9:25S-02  1.84E-01  4.56E-01  1.63E-00  5.2JE-00"li25E+01 
2040.  6.31E-03  1.63E-0B  1. 577-05  2.54E-03  6.41E-03  3.08E-02- 8.30E-02  . 
2090.  7.4PF-A2  1 • 1SE-C1  2.56E-01  5.09E-01  2.06E-00  6.10E-00  1.37E+01 

2100.  1.81E-03  l'.47E-0t  2.10E-03  2.65E-03  1.06E-02  4.28E-02  1.08E-01 

2100.  7.89E-02  1 .675-6 1  4.95E-41  6.00E-01  2.54E-00  7.00E-00  1.47E+01 

2110.  9.55E-04-  1.26F.-03  2J32E-0J  1  .'S9E-03  U73E-02  5.86E+02  1.41E-01 


2110.  8.13--02  2.01E-0!  5.00E-01  7.60C-01  3.04E-00  7.90E-00  1.576+01, 

21201  5.652-04  S.B2E-04  I.70E-03  2.46E-03'  2.80E-02  8.0SE-02  1.81E-01 

21201  6.4OE-02  2.14E-01  5.58E-01  9.73F-01  3.5SE-00  8i76E-dO  ,1.6'6C+01 

?J3Qi.  4,882-0.4  6.'96E-Q4  1.43E-03  4.98E-03  4.38C-02  1115E-01  2.31Ef01 

2130.-  6.611E-02  Zi35£-01  6.12L-01  1.24E.-C0  4.19E-0Q  94602200  1.75E+0.1" 

214&;  4.78E-04  7J30E-04.  1.  73Et03  S.S5E-03  6.73E-02  1.63E-0J-.-2.95E-0.1 — 

“214Q.  Si97E-02  2. 508-01  6.6.6E  +  01  l.frOE+OO  4.85E-D0  1.04E+01  1JH2E+01 

2160.  4.91E-04  7. 8 1 £-04  2.20E-03  1.326-02  1.016-01  2. 256+61  3.7JE-01 

2150.  '}.2aC-02  2.66F-01  7.40ErOi  2i0iE-00 -5.'63E-00;?l'i'12E+Ol  1.90E+0T 

2160.  5.22E-04  8.566-Q4  r.7?c-01  Bi'4.6+02  'ii$3£-01  3.06E-01  4V68E-01 

2169.  9 ,'605-02  21456-01  «.'77E-0l  2 .516-00.  6.486-00  1.206+91  .U986+01 

2170.  i'.ait-Qi-  +.O?E-03  t'il6E-02  2.2.8E-01  ^.{lTErOl  5.80E+01 

2170.  9.,37C+02  liBSC+O.l  1 .  !  1E-0O  35 13E-00  -7.J5E-0C  1 .3CE401  2.032+01 

.ISO.  6.9Si.-0«1.46E-0'3  8.8?Er03-  K07E-d'  3V32E-01  5436E+01  7*106+01. 

2-1  BOi-  il. 0*6-01  4.22F+ai^l.4oE-0Q  »; 9ie-00  8126E+Q0  .li39E+Ol.-2i05E+Ol. 

2  TOG.  OiOlErOA  2.60FrO3  2+18E-02  1  iS&e^Ol  'R70E-dl  7i0>E-01  ’8J57E.-01. 

2190.  *.l?T~il  4i90E-01  !»72E“W4,8lE-00  °»27E+00  1.4.7E<01  2.06E+01 
2260.  1.2H6-03  6i27t-03  .5, 0CE-02  3.152,-01  6. 56E+01,  3;o3E-01  l.OOE-OO  • 
2.90.  1.296+0,1  5.816-01  2.1,0‘E-<>p  1.81C;ftp  4-i05t+01  .li56F.+0l' 2.1PE+01 

2210.  2.45C-03  1  .65E+C2  i  .OOl-Ol  4;92E"01'  8i88E-01-  1.11E-00  1.176-00 

P21 t.  1.525+01  6.97E+01  215'E-OO  ftiOCE+GO  U19E+01  . 1 .64£+ol  2. 17E+01 

2220.  7 .,9 >£-03  4.54S-08  P.OiE+fW  7;28'E+01  1.166-00  T;35E+00  1.36E-00 

2320.  1.577-0’  6i4i*F+ol  liC2C-00  8.028:00  1.33E+.01,  .1., 726+0?.  2.25E+01 

2230.  3.10E-02  1.22E+01  J. 672-31.  lVC'ES-eo  :1.4  6E-7)0  i.61E+00  1455E+00 

2230.  2.65E-A7..  M-iHE-C'C  3;6»E-00  9i  12E-d  li,45v<0!  1. 8.16*01  2.35E+01 

2240;.  '*.845-0*  /.66E-31,  7.05t+ni  1  .466-^0  -U'8»cr,00'  4,86Et/M  1.72E-06 

2;40l  -.J'ErOl  liieiF-OO  +.445-o>  lidjC+oj  r.JisF+Oi  l:.89f+()3..2.41E+01 

2250.  2.32E+1-1  -i«6E-01  5.i3E-r.!  !.9SE-00  -2 i22E-Of*  ,2.KE4oQ  li?7E+0C 

2250.-  6.-t»3£-o:  a.y'iEf'tO  O.’Pit+OO  l.ME>oi>  :r.5?c+Cl’  1  SjB^C+6yl‘.2i 3^E+pl’ 

2260.  H-m-Qi  7im-05  i;53t-iOP  ,c;65£-dC' 2.tlE-90  ;2i‘3'2Er,00  li9SE?:00’ 

2260.,  1.015-OQ  2i6^t-6o  5i99E4>),0  l’ilOE+Ol.  i'.t9E+01  l'.85E+Ol: ^ilJt'+O.l’ 

2270.  5.2rEr.01  1.26E-00  2ij4Er0t  3i27E-Oo  2^96c-.C0  2i*9E+bO  2i06E?50  . 

aPO..  U  jaS-rOO  3.2CE-0C  .  i55F>00  ■-i.26c-+aj  '  1.5.7E+61  l.'78E+b:.-'U96E+01 

33Sbi.  7.90Erftl  Zi^E-O'i  3V525-00  3,947+oi.  .3^5^00  2i6iF+i);)  .2.i.0E-00 

2280.,  liill— 0'  3376E+0.4  J34,90E-h+  U27F+01  ftC52F,+bl  .loTEtOlt  1474E+01' 

224oi  Ii84j-0O  -S'.lOCrOQ  4.53E-00  V.48E-00-.2'.‘69E-00  .2.12E:00. 

_?•.?£<?*•  2iS/t-njl’«Jl27E-0i}  7.hOt-DO  r.232’-Ql  i»4BE+0r  1.58E+G1  l.'59,E+01. 

iioO.'  r-4-.4'OE~00  :7.42r-0<i  67v7?-0b-  V.VSE-OO  3.60^0  2.71E>00  2.!lE:0C. 

_  ?4<jn.  iiOlF-OO  5.385-400, fliVorr?'/  li27E+0U.i.42E+Cl  :i.47E+OM>45E+01 
2316.  6.72>--(10  1. 1 2F+I7S  7.905+00  i.'14r4no  3.631-69  2.6OE-.90  2407E-00 

23roi  SiTTC-E-AO  7 J395*bC.Ti0ib+Oi-^2JE3 01' +01-  •l".!3*E*b l- ■  1736E+01- 
2320.  liJ4ctrf>  1 .465-07  S.P5C-00  6.l4Erno  3.60E-CC  2.65E-00  2.01C-00 

r< 26.  6.j>E-:nf  ItO.OE+Ol  l.tZUal  » i  1  feL  +A 1  1.226+01  1.21E+01  lilsE+Ol 

2330.  1.6/c+Ol  1.4.:F»Cl  Z.97E-00  5.08E-00  3.57E-00  2.61E-00  3.95E-00 

23?0.  li2.'++01  1.145-40;  r.l2E+'0V‘ 1.08E+0'.  li'10E+91  1.07E+01  l.OOE+Ol 

2340.,  l.'.ToC *91  l.UE+01  7+4'£e4qi)  6,U,E-:r.o  3.55E-00  2.63E-00  l,;25E-00 

23,4(1..  J.Ate+oi  l.,105+*M  r.Ctfc+01  8.E5E-OC  ^i70:+0C  9.2-E-OO  8.52E+00 

2944.  i.73r+ftl  ^1-j-tOl  7.6(7+00  5i27E-0(>  -3F50C-*b4  ZJ40E+00  i.S9E-00^ 

2;.30.%  1.64F,-9l  t.'.n7',ioi  ;.i6E-tl0  8.51E-(iO-  “VlbEr0C  7i61E-00  7;7.2E-0'' 

2S*«" 1.70F+O1  1.66*-i‘01  1.00.-‘01  5.3Jc^do  3.'28E-C0  2.126-0,0  1-.43E-00 

2366.,  l.ant+ol  C.^IE+OO  7.606-00  6.5TE-0O -6.33E-00  ;5»7Y5-00  5.45E-00 

2-376.  2..  E+ri  ..‘a0E*01  /.50E-»r.  4.9fe£TC0  -2i.69EfO0  1.67E-0O  1409E+00 

237.1.,  fr.v.7t+t.O  +.4 .2-00  5.V0£~y  a.tOt-OO-^iibE-OO  4.12E-00  ^.67E-00 

PJcO.  liOiEKii  7.2oF-«i-u..vo£-0m  3.57;--  <'C  Ii77E+00  Ii05c-0C  6150E+0I, 
21«0.  2.63E4.-.0-  3Y;-.'-CC  4.10E-'O  2.8CE-00  2.PCE-O0  1.90E-00  1.60E-00 

+30,,.  2;-'7t  -cl  4.+7^-Ox  1.70f->>->  1.207-00  7. FOE-01  4.P0C-O13.O0E-61 

2  +  4...  7.  +  .  2-01  1.1-7+.  l.ofig.nn  1. 49;2O0  1.20E-70  V.0.96— Of.  9.00E-01. 

2470.  6.y>E2f,2  *.«  >.*  -uj  l.flOE-n;  1.506-9.1  U50E-01  1.30r-0~l  l.Ob't-Ol 

;-4jrt.  4,y</3pl  7..i:-Ji  6.63^-01  (.062-01  S.60£-7l  87'OE-Ot  8iO0E-Ol 

-.faMArii:.  F.A  BAOIAI Itin  2A0IAII0M  KAOlAIltlli  RAOiAIlON 

2A1IAI10'  R7.4|A1‘CS  kAO  I  A?  lo\ 

ft  MAH  70 

n3|<i.  >U»ft01ET  Alr.C.lAfT  -  iE  LEVEL  COimHIONS  -  900  PFO  EC>T  STATIC 

7  ■ 
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are  finished  for  a  particular  ALPHA,  a  (2)  will  return  the  program  to  accept  a 
new  ALPHA  and  ISTOP,  and  a  (3)  returns  the  program  to  read  a  new  set  of 
flow  field  data.  The  input  data  sequence  is  summarized  in  Table  V. 


TABLE  V.  INPUT  DATA  SEQUENCE 


Card  No. 

Column  No. 

Description 

Format 

1-72 

1-7 

*  Wavenumber 

F7.1 

1-72/odd 

8-70 

* 7-S/d  values  (300,600 . 2100eK) 

7E9.2 

1-72/even 

8-70 

* 7-S^2/ d  values  (300,600 . 2100 *K) 

7E9.  2 

73-74 

1-72 

*  Printout  headings  (RADIATION) 

12A6/4A6 

75 

1-12 

Date 

2A6 

76 

1-72 

Title 

12A6  ' 

77 

1-5 

Number  X  Stations 

15 

6-10 

Number  Z  points  each  X 

15 

78 

1-10 

Aspect  Angle 

F10.5 

11-15 

ISTOP 

15 

*  Standard  for  all  runs. 


A  typical  set  of  output  data  from  the  program  is  shown  in  Table  VI.  The 
aspect  angle  is  listed  in  degrees,  the  downstream  position  is  given  in  centimeters, 
and  the  plume  thickness  at  the  downstream  position  (X)  is  also  given  in  centi¬ 
meters.  At  each  downstream  position,  the  thickness  along  line  of  sight  is 
divided  in  a  number  of  segments  to  make  the  radiation  calculations.  The  number 
of  such  segments  or  zones  is  denoted  by  STEPS  and  the  size  of  each  step  is 
indicated  by  DELTAL  in  centimeters.  These  quantities  are  automatically  cal¬ 
culated  by  the  computer.  The  Z  coordinate,  previously  described,  gives  the 
radial  position  at  which  the  calculation  is  being  made  at  a  particular  downstream 
location.  This  RADIAL  POSITION  is  also  given  in  centimeters. 

The  first  two  tabulated  columns  of  data  give  the  wave  number  (cm-1)  and 
wavelength  (p).  Next  N  radiation  columns  are  listed,  of  which  the  first  N-l  gives 
the  radiation  emitted  from  the  plume  surface  in  the  direction  of  at  the  down¬ 
stream  location  (X)  and  radial  position  (Z).  This  is  given  in  watts  per  steradian- 
micrometer-cm2.  The  last  radiation  column  for  which  no  radial  position  is  indi¬ 
cated  gives  the  spectral  radiation  per  unit  iength  of  plume  surface,  i.e. ,  the  first  N-l 
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NOTE  -THE  LAST  AADtATlO*  COLUAN  AAOVfc  GIVES  TH€  SAECtAAl  WAPI A f T  LtHGlH  Of  fLUmf  SURFACE 


TABLE  VI.  TYPICAL  OUTPUT  DATA  (Continued) 


'  ~  **  ~  5  *Wi  R*TOP!^5!55^R^TO?5^5?!^; 


3®?5?5? 


values  of  radiation  have  been  integrated  across  the  radial  position  coordinate. 

The  last  column  therefore  has  units  of  watts  per  steradian-micrometer-cm.  The 
next  to  the  last  line  of  printout  gives  the  radiation  per  unit  length  of  plume  ■ 
surface;  i.  e. ,  the  last  column  has  been  integrated  over  the  spectral  regime. 

As  indicated,  this  has  units  of  watts  per  steradian-cm. 

After  the  spectral  and  spatial  distributions  have  been  computed  for  each 
downstream  location  at  a  particular  aspect  angle,  the  results  are  then  integrated 
over  the  whole  surface  of  the  plume.  The  last  output  page  for  a  given  aspect 
angle  contains  this  information  as  a  function  of  both  wave  number  and  wavelength 
This  is  listed  as  SPECTRAL  RADIATION  and  given  in  watts  per  steradian- 
micrometer.  Next,  a  spectral  integration  is  performed  which  gives  the  total  radiation 
emitted  in  watts  per  steradian.  Finally  the  centroid  of  radiation  along  the  length 
of  the  plume  is  computed  and  listed  in  centimeters. 
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Appendix 

FLOW  FIELD  PROGRAM 

This  appendix  outlines  the  flow  field  frozen  mixing  computer  program 
which  is  used  to  provide  the  input  data  to  the  radiation  program.  A  typical 
set  of  input-output  data  is  shown,  a  description  of  the  input  data  is  given,  and 
the  computer  program  is  listed.  A  brief  description  of  the  viscosity  options  is 
also  included  for  completeness. 

The  viscosity  option  occurs  on  the  second  data  card  and  is  labelled 
ITURB  by  the  computer,  where  ITURB  =  0,  1,  2,  3,  or  4.  The  statement 
ITURB  =  0  indicates  laminar  flow  and  the  viscosity  is  given  by  Sutherland's  law. 
For  ITURB  =  1  or  2,  the  eddy  viscosity  is  computed  from  the  expressions 

H  =  K2  r,.  p„Uo  ,  (ITURB  =1) 

•1 

M  =  K2  r^ jp0  Uo  -  Um  |  ,  (ITURB  =  2) 

where  K2,  the  turbulent  viscosity  coefficient,  is  taken  to  be  0. 0285.  The 
subscripts  refer  to  jet  center  line  and  ambient  conditions,  and  rly  is  the  value 

I  1  I 

of  RatpU  =  -Pouo)-  For  ITURB  =  3,  r,^=  |  r0>99  -  r0<50|,  where 

r„  50  is  the  value  of  R  at  U  =  1/2(Uoo  +  U0^  and  r0  99  is  the  value  of  R  at  U  = 

(o.  OIU^  +  0. 99U0).  The  viscosity  is  then  given  by  the  expression 

H  =  K2  rv  Po  U„  .  (ITURB  =  3) 

'2 

For  the  option  ITURB  =  4,  the  viscosity  is  computed  from  the  expression 

n  =  10"4  +  X(pw  Uw  +  P„  U„),  (ITURB  =  4) 

where  X  is  the  downstream  position.  The  latter  option  is  utilized  for  jet 
aircraft  plumes  in  the  core  region  with  this  region  being  defined  by  the  condition 
d2U/ dR2  =  0  .  Outside  the  core  model  3  is  used  for  aircraft  in  flight  and 

R  =  0 

model  1  for  static  conditions.  The  program  automatically  selects  the  mode 
outlined  above  when  the  input  data  specifies  option  4.  Table  A-I  describes  the 
input  data  sequence  for  the  program.  The  thermodynamic  data  are  obtained 
from  NASA  document  SP-3001,  pages  308-326.  A  typical  set  of  input  data  is 
shown  in  Table  A-II  followed  by  the  program  itself  in  Table  A-III.  Table  A-IV 
then  shows  a  few  pages  of  selected  printout  or  output  generated  by  the  input  data 
described  in  Table  A-I.  The  initial  page  specifies  the  input  data;  the  remaining 
output  gives  the  properties  as  a  function  of  downstream  and  radial  position. 
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The  onlv  programmed  error  message  occurs  if  the  summation  of  the  initial 
specie  mass  fractions  (indicated  by  SIGMA)  differ  by  more  than  1  percent  from 
unity. 


The  actual  output  from  this  program  that  is  used  in  the  radiation  model 
is  written  on  magnetic  tape.  The  program  first  expands  the  number  of  radial 
points  to  25  by  defining  the  plume  width  to  be  that  radial  point  where  T/T^  =  1.05. 

This  width  is  then  divided  into  25  evenly  spaced  points  and  the  properties  deter¬ 
mined  at-eaeh  point  by  interpolation.  Dimensions  are  converted  from  feet  to 
centimeters,  mass  fraction  to  mole  fraction  or  partial  pressure,  and  temperature 
ratio  to  absolute  temperature.  The  resulting  tape  is  then  read  as  input  to  the 
radiation  model. 


TABLE  A -I.  INPUT  DATA  SEQUENCE 


Ca  rd 
No. 

Column 

No. 

Description 

Format! 

1 

1-72 

Title  Card 

72H 

2 

1-5 

Month  (Jan  =  01, . Dec  =  12) 

15 

(5-10 

Day 

15 

11-15 

Year  (Last  Two  Digits) 

15 

1(5-20 

Initial  number  of  grid  points  (if 
the  variable  profile  option  is 
used)  or  initial  number  of  points 
in  jet  (if  the  step  input  option 
is  used) . 

15 

21-25 

Number  of  species  (21  maximum) 

15 

26 -JO 

Pressure  option: 

0  =  Constant  Pressure 

1  =  Polynomial  Fit 

15 

31-35 

Viscosity  option: 

0  =  Laminar  (Sutherland's  Law) 

1  =  K2r^  p0U0 

2  =  K2r«/2|^Uo  *  ^eUe| 

=  K2ry  p0Uo 

■l  =  10-1  +  X(p0U0  -  Peue) 

15 

:i(>-  to 

Flow  option: 

0  =  Axisymmetric 

1  -  Two  Dimensional 

15 

TABLE  A-I.  INPUT  DATA  SEQUENCE  (Continued) 


Card 

Column 

No. 

No. 

Description 

Format 

41-45 

Input  profile  option: 

15 

0 

=  Step  Input 

1 

=  Variable  Profile 

3 

First  print  increment  (ft) 

Final  X  for  first  print  increment  (ft) 

Second  print  increment  (ft) 

E10.  8 

iH 

Final  X  for  second  print  increment  (ft) 

E10.8 

41-50 

Third  print  increment  (ft) 

E10.8 

51-60 

Final  X  for  third  print  increment  (ft)  and 

for  terminating  the  case 

E10.8 

61-70 

X  initial  (ft) 

E10.8 

4 

1-10 

Lewis  Number 

E10.8 

11-20 

Prandtl  Number 

E10.  8 

21-30 

The  initial  radius  of  the  jet  in  feet  if  the 

E10.  8 

step  input  option  is  used,  otherwise 

A  the  radial  spacing  in  streamline 
coordinates 

31-40 

Kj  the  coefficient  of  laminar  viscosity 

E10.8 

41-50 

K2  the  coefficient  of  turbulent  viscosity 

E10.8 

5 

po 

Pressure  fit  coefficients: 

pi 

If  option  1  is  used,  P0 

E10.8 

is  the  pressure  in  Ib/ft2 

21-30 

P2 

and  P^  P2,  P3,  P4  are  blank. 

If  option  2  is  used,  Pj  is  the 

E10.8 

31-40 

P3 

coefficient  in  the  pressure 
polynomial  as  follows: 

E10.8 

41-50 

P  (lb/ft2)  =  £  p  XI 
i=  0 

E10.8 

6 

■ 

Thermodynamic  data  for  the  different 
species  being  considered.  There  are 

3  cards  for  each  specie  as  described 
below.  The  specie  may  be  input  in 
any  arbitrary  order. 

6a 


f 


1-6 


Hollerith  representation  of  the 
specie  (e.g.  H20) 


A6 
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TABLE  A-I.  INPUT  DATA  SEQUENCE  (Continued) 


Card 

No. 

Column 

No. 

Description 

Format 

mm 

Molecular  weight 

nm 

B 

T 

1L 

lower  temperature  bound  for 

fit  1 

3 1  -40 

T 

1H 

upper  temperature  bound  for 

E10.8 

fit  1 

41-50 

T  l  lower  temperature  bound  for 

E10.  8 

fit  2 

51-60 

T  upper  temperature  bound  for 

E10.8 

fit  2 

6b 

1-10 

ai 

E10.  8 

Coefficients  for  low 

11-20 

a2 

temperature  fits 

E10.8 

21-30 

a3 

E10.  8 

See  NASA  SP3001 

31-40 

a4 

E10.8 

41-50 

a5 

E10.8 

51-60 

a6 

E10.8 

61-70 

a7 

E10.8 

6c 

1-10 

•• 

Coefficients  for  high  temperature 

E10.8 

11-20 

fits 

E10.  8 

21-30 

E10.8 

31-40 

NASA  SP3001 

E10.8 

41-50 

E10.8 

51-60 

E10.8 

61-70 

E10.  8 

IF  THE  VARIABLE  INPUT  OPTION  IS  USED: 

7 

1-10 

To 

Axis  value  of  temperature  (°K) 

11-20 

T, 

E10.  8 

11-20 

• 

The  values  of  the  temperature 

E10.8 

11-20 

. 

at  each  of  the  inputted  psi 

E10.  8 

11-20 

. 

grid  points  are  punched 

E10.  8 

61-70 

T  6 

7  to  a  card  from  the  axial  value 

E10.  8 

and  ending  with  the  free  stream 

values  (°K) 

7' 

t7 

E10.8 

t7 

1 

E10.8 
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TABLE  A-I.  INPUT  DATA  SEQUENCE  (Continued) 


Card 

No. 

Column 

No. 

1 

61-70 

1-10 

00  00  00 

_ I 

1-10 

1-10 

61-70 

1-10 

1-10 

1-10 

61-70 

1-10 

9 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

1 

Description 

Format 

Tv 

E10.  8 

Tl3 

Tl4 

E10.8 

Tl4 

E10.  8 

Th 

E10.  8 

T 

Free  stream  value  of  temperature 

6 

(°K) 

u„ 

tmm 

Uo 

The  values  of  the  velocity  at 

E10.8 

each  of  the  inputted  psi  grid 

E10.  8 

points  are  punched 

u6 

7  to  a  card  from  axial  value 

E10.  8 

and  ending  with  the  free  stream 

Uv 

values  (ft/ sec) 

E10.  8 

E10.  8 

E10.  8 

U,8 

Ul4 

E10.8 

Ul4 

E10.8 

Ul4 

E10.8 

U 

e 

Free  stream  value  of  velocity 

IF  THE  STEP  INPUT  OPTION  IS  USED: 

ai 

«2 

illiliS 

«3 

Axis  values  of  species 

E10.8 

«4 

mass  fractions. 

E10.  8 

E10.  8 

E10.8 

«7 

E10.  8 

«1 

®2 

Values  of  species  mass  fractions 

E10.  8 

«8 

at  each  inputted  psi  grid  point 

E10.8 

«  4 

are  punched  on  one  card  per  point 

E10.  8 

<*5 

beginning  with  the  axial  value  and 

E10.8 

«« 

ending  with  the  free  stream.  A 

E10.  8 

«7 

maximum  of  7  species  per  card. 

E10.  8 

a7 

The  order  of  the  species  must 

45 


TAH1.E  A-I.  INPUT  DATA  SEQUENCE  (Concluded) 


Card 

No. 

Column 

No. 

Description 

Format 

(51-70 

o7 

correspond  to  the  order  of  the 

(51-70 

°'7 

thermodynamic  data  (card  type  6) 

<>" 

1-10 

"I 

BED 

11-20 

O' o 

1911 

2i-:io 

«'3 

E10.8 

:u-40 

°'4 

Free  stream  values 

E10.8 

41-50 

"5 

E10.  8 

51-60 

«'6 

E10.8 

(51-70 

°'7 

E10.8 

IF  THE  STEP  INPUT  OPTION  IS  USED: 

7 

1-10 

Jet  temperature  (°K) 

EH 

11-20 

Free  stream  temperature  (°K) 

■ 

Jet  velocity  (ft/ sec) 

E10.  8 

■BEK 

Free  stream  velocity  (ft/sec) 

E10.8 

8 

1-10 

°t 

Specie  mass  fractions  in  the 

mrwm 

11-20 

O' 2 

jet.  A  maximum  of  seven  species 

■  - 

«3 

per  card.  If  there  are  more 

E10.8 

.*11-40 

than  seven  species  continue  on  another 

E10.8 

41-50 

«5 

card  in  the  same  format.  The 

E10.8 

51-60 

order  of  the  species  must  corres- 

E10.8 

61-70 

O’ 7 

pond  to  the  order  the  thermo- 

E10.8 

dynamic  data  (card  type  6)  for 

each  specie  input. 

5) 

m 

O’l 

Specie  mass  fractions  at  the 

«2 

edge.  Comments  for  card  8  are 

■iiilll 

■ns 

O' 3 

applicable. 

E10.8 

.‘11-40 

°’l 

E10.8 

11-50 

0  5 

E10.8 

51-60 

O’fi 

E10.8 

61-70 

O' 7 

E10.8 

10 

mem 

The  initial  radius  of  the  jet  (ft) 

E10.8 

■89 

Run  number 

15 
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ABLE  A— III.  FLOW  FIELD  COMPUTER  PROGRAM 


MAlNl  -  FfN  SCONCE  STATEMENT  -  1FMS1 


U  3/09/70 


C  HA  IN  OfafaStOlO 

C)m«I:N  R.X.OPOX.POOT  .HI  ,FE  ,DFUUl.F,l)lLPSI  .  CX  ,XMUl .  XMUL  .DEL  .XML'K  l.X 
lwtlR/  .PkM  .PCM  ,UF  •  RHCE  .TE.XP3.RFALF 

c  jur-t  itiut.HPsr.MHAir.KUiT.  if  i  ms.kpsi.i  page,  i  press.  I  start,  ims 

1.1  lURA.IJEt, NS, Ml>AY. MONTH, MYEAR.1AX2C.IVAR 

c'im-cn  NTMoi.ei?i».c(iEffp(s i,AiPHAia,59i,HcuTr5vr,'Kcarms>r.irrr5<>r'  ■  " 

l.T  I  MSI  ,M(  21  .59),  RALPHS  <2  l.59I.PS  I  1591.  PUT  69).  UC  59) , XLEl 59) , t| TLEI l 
22)  .SUM(S4).FTA|59)  .XM0159)  .Also)  ,A*AI  M59)  .RHCOUT ( 59 ) , UCUTI59)  ,7()U 
)M  E9  I.CPi  2l.5<>|,HSIG(59)«AHCi59).Yt59),XMAXI)) ,  X  FHT  (  )  )  ■  CP8 AHI  59 )  ,  S 
A  tr.MAlSSI.HrtAR  <«<>).  AH  T  SI  21  ,2.7)  .tf  1TS( 2 1 , A ) . SPEC  101 21) 

HFttlM)  10 

1  Call  input  . 

C.ML  IMTAl  JANS  10  33 

2  CALI  FKOM 
CALL  CCMPUT 

it  hums. i.t. i)  ou  in  i 

CALI  fcXPLtC 

1M  IHMS.GI.il  GO  TO  1  " 

GO  IT  2  066S10A7 

EMC  .  "  066S10A8 

03/00/70 

FRO/?  -  FFN  SCLRC  E  STATEMENT  -  IFMS1  -  '  '  "  " 


1 

3 

5 

7 

9 

l* 


XUt'P  Gill  INF  INPLT  ObbSOABO 

Cu  4MCN  K.X.OPCX.POUT.HE.PE.CPOUT.P.OELPSl  .CX.XMOT, XMUL. DEL  .XMUK1.X 
IMu^.PRNT.PCNT.UF.RHCE.TE.XPJ.RMLF 

CUYfCN  lOUT.HPSI .MHALF .MINI T. IF  I M S. KPS  I . 1  PAGE . I  PRESS. I  START, tTURft 
l.l  IlJRA.  1JE  T  ,  NS.  MOA  Y.  MONTH.  MYEAR  ,  I  AX2C.  I  VAR 

Common  RTMulE  (?1 1  .COEFf  PIS  ),  ALPHA  l21.59),HCUTU9nTC0TTmrnm39) . 

l.T  (5  9)  .HI  21  .5S),RALPHA(21.59),PSI(59>.RUt59).UI69),XLE(59),TITLEt  1 
2  2)  .SUMlS9I.FTAt55),XHU(59l  , AI59  ) , ARAI T 159)  .RHG0UTC59)  .UUUTI59 )  VTOU 
)T(  ‘4  ).CP(  21 .59I.HSTr.(59I.KHC«S9»  ,Y{59)  ,XMAX  1 3) ,XPRT  1 3)  .CPbARl  59 1, S 
AK.»<AI5S).HBAR154I.AFHSI21,2.7)  .  !F  I TSI 21 ,4 )  .SPEC  ICC  21)  ' 

01‘FNS  1«1N  YSJET  (2  1)  .YSEUGE  I  21  I 

IF  IN  IS *0  - - 

IP  Ai,E*G 
It)  Li  *0 
CX  .U  .0 

KO,  57529  32  '  066S0502 


RfA, MS. 1)  I  T  ITLEI  1 1  .1*1,12)  3 

RFAOIS.IOOImCNTh.MOAY.MYEAR.MPSI.AS.IPRESS  ,  I  TOR  8  .TAV  20  iTVAR  . .  *—  12 

NPSI*MPSI-1  Q66SC528 

Rf  ADlS.lOiA  IXPRTIl  I.XMAXU  ),XPRT  12  ) ,  AMAX  1 2  )  ,  XPRT  Cil.XHAXOt.X  22 

RCAiMS.lOCCI  XLEII).S!GMA(1),YJE1.>MLK1,XMLK2  23 

RCAOIS.IGCO)  ICtlEFFPU  1 .1*1.5)  2A 

00  SO  J*l.NS 

RF.MMS.200C)  SPECIOI J)  .  WTMOLE  I J  I  ,  I  TF  ITST  J  ,K)',K*1 ’.AT)  33 

RF  At)  1 S.  1000  )  (AH  ITS!  J.  l.KI  ,R*1  .  7)  AO 

RCAOIS.10C0)  1AF1TSI  J.2.KI  .Ml. 7)  S5 

*0  CuM  INUE 

1FI1VAR.NE.C)  GO  TC  AO 

RCAtMS.lOOu)  T  JE  T.  TEOGE  .  C  JE  T  «UE  CCE  55 

RCA0C5.1000)  lYSJET(J)  ,J*l.NS>  56 

re  ao 1 5 • looo )  ( y sedge ij  i  , j* i *nsi  63 


CUM*C.C 

P*COEFFPI  1 )  ' 

Co  )  7  J*1,NS 

37  CUN*CUM»YSJtTU)/rtTMOLElJI 

RHUJ  .P/84517. 5C1/TJFT/0LM  *  . 

IF  « I  AX20.C0.U  )UPSl*YJET»SGRTIRMCJ»UJH)/FlCAT»MPSt-l)  79 

IF  I  I  AX  20.ro. 1  I  Ocl P  S I  * Y JC  T*RMO J  *L  JFT/FLGAT IMPSl“l J 

IX)  3P  I  *  1 . M P S I 
Tt  I ) *T  JFT 
U  1 1 *U JFT 
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TABLE  A— III.  FLOW  FIELD  COMPUTER  PROGRAM  (Continued) 


- ro-Tr'jii;Ns-  —  - - 

36  AL PHA!  J» I ) -YSJETIJ  > 

MPSI  -MPSI  *2  - - 

NPSI-MPSI-1  Q60S0114 

TIMPSIl-TEOGE  - -  -- 

TINPSI  l-TEOGE 

- urm>ni<OEOGE  ■  •  : . -  •  - - - 

UINPSI  l-UEDGE 

CO  39  J-l.NS  - 

«l PHA! J.MPS II -YSEUGE ( J I 

~39'ALPHA(  JtNPS  I  l-VSEOGEi J  )  - 

GO  TC  60 

— *0  DEEP  ST-YJET  - - 

NPSl-MPSI-1 


03/09/70  | 

FRO  12  -  EFN  SOLACE  STATEMENT  -  ‘TFKT5T- -  '  ~  '  | 

- ifFitnynocoTTTur.tii.Mpsn - tit 

RE  AD  (St  1000 )  «U(  l>. 1*1, MPSI)  121 

CO  41  I*  1*  MPSI  '  . . .  ’  ' 

REAO(StlOOO)  ( ALPHA!  J,  I )  ,J*1  ,NS )  13C 

'41  CONTINUE  . . 

60  RETURN 

“1”  FDKfflT  1 12A6 )  .  :  oSnS0646 

100  FORMAT ( 14151 

1000  FORMAT  ( 7E 10 .8 )  - OWTSObW - 

2000  FORMAT!  A6t4Xt5ElC«  8) 

"  ENC  - -  06650649" - 


FR023  -  EFN  SOURCE  STATEMENT  - — fPKT57 


SUBROUTINE  FRCZ1 

COMMON  Rt Xt OPCXt POUT .HE,  PE  .OPOUT.P.OELPSITCX.XMUTTxAULTUEI.- .XMUKT.X  ” 

1MUK2  tPRNT  tPCNT  tUEt  RHOE  t  TEt  XP3t  RFALF 

COPMCN  10UT.MPSI  .MHALF  .MINI  T.1F  INI  StNPSrn^ETl'PSISSTrSr/rRTTrTTrRB - 

It  1  TURAt  IJET .NS.MOA V.MOMH.  PVEAR  . I AX2C.  IVAR _  _ 

1.T159)  tH(  21 . 391  .HALPHA  (21,59),PS1(59),RU(59),U(  59)  t  XLE  (59), TITLE!  I 

22)  .SUM!  59),  ETA!  59)  ,XMU !  59)  .  A 159 ) ,  A  MAI  TtS^nKBUgUT  1 59)  .uOuT  139).  TOu - 

3TI 59 ),CPI 21 .55) .HSTGI59 ),RHC(59),Y(59) .XMAX (3) . XPRT 1 3) .CPBAR! 59 ), S 

4lGMA(55).HBAR(59).AFITS(2l.2.7),Tfirs<  2T74T  .SPICT!)  (Tf) - 

P*CUEFFP! 1 ) ♦X*(COEFFP! 2  )+X*(CUE  FFP (3 )*X* (CCEFFP! 4)+X*CCEFFP(  5 ) ) ) )  066S1396 

DmiC0EFFP12)4X*r2.6C0EFFP(31>T*nV«CPEFFPI4)»X*C0EFFPI5)64.)) - D6A51397 

00  1  I  *1. MPSI 

ROCTT!  I)*SOKT(T(  1) )  - - 

AM  AIT!  I )  *  0 .  C  066S1440 

CO  2  J*1,NS  - - 

2  AW  AIT!  I)*AmAIT!I  )♦  ALPHA!  J,  D/MTPCLEIJ) 

'I  RHOI  T)«P/89517.5C1/T!I )  /ANA IT ITT - 

DO  20  I*  1, MPSI 

DO' 1C  J-l.NS  - - 

K«1 

IF(T!I).GT.TF1TS(J,2))  K»2  . . . . 

CP!J.  I  )*(  AFITSIJ.K.  X )  ♦  T  <  I » *  ( AF  ITS!  J.K.  2)*TU  )*! AFITSIJ.K,  3)*T!  I  )*( 

lAFITSf  J,A,4)*TII  )*AF1TS(  J.K.5) ) ) )  )’*l79fi7'i'l  .S/wTMCLE  <  J  J - 

HI  J.  I) -(AFITSI  J.K,  6)4T(!)*!AF  ITS!  J.X, l  )»T  (.!)*!  AFITSU.K. 2) /2.4T1  I) 

l*(AFiTSIJ,K,3)/3.+  T(I)*!AFITS!J,K,4)/4,+T(l  )*A’FITS(  J,"K, 5T/5T) ) )’) ) - 

2*1  .9  87*  1.8/mTMOLF!  J) 

10  CONTINUE  "  - - - 

CP  OAR!  U-0.0 

—  "HBJRin-0.0  - 

CO  3C  J-l.NS 

CPUARl  II-CPBAR(!)*CP(J,1)*AIPHAU,I)  .  ’  - 

30  H3ARII  l-HBAR!  !)»H!  J,I)*ALPHA!J,  I) 

20  CONTINUE  '  - - -  ----- 


Yl  11*0.0 
WT2?  1-2, MPSI 
IF  ( I  AX  2O.EO.0) 


066S1535 


'CA0S02A0 


49 


I  ABLK  A— III-  FLOW  FIELD  COMPUTER  PROGRAM  (Continued) 


I V i  !  I  >SOMTIV(  IPSI»IPSI  I  I l/KMflll  l/Ul  IIaPSII  t-l|/RHOIl-ll/U060S02<.l 

2tl-tl!l  88 

'*  I  I  AX  2u.  *-0  .  I 1 

iri  II  *YI  I  -  1 1  aCUPsI  *1  l/ul  1  lAl./MMUl  l-ll/UI  1-tl  1/2. _ 


28  CoM  INUF 

l>  (  I  tiJKrt - f.e-  .Cl  uC  I"  4C 

xrM  •  1.0 
cn  23  im.rpsi 

? i  xnki  i i-xmuk  i*iii  i •  i»t;ci mi /« i( i  i«m.i*xniT  066SIS49 

X«U  *'.MU|  ti  _ _ 

Gil  fO  tco  "  "  666Sl(J^6 

4C  Uo  101?*..  <.8,6C,IJ01I.!  TU«A 

<.8  lu  M*  .  8  •  I  w  miM  1 I  Ay  1 1  I  ♦Hull  IMPS!  )*U(PFSI  1 1  066S1558 

L!  82  J*l.MPjl  066S1559 

t » Ml'  SI -J»  t  066S156Q 

If  I-  hill  I  »«Ul  t  1-  )UM  I  81,81,82  066S 1561 

8?  ConTInuF  *  046S156T 


91  CU  *Yl  l  I-  (  V  I  I  l-Yll  ♦ll)»(HMCUI*llll-EUH)/»#nCUI*UI  ll-RMCl  l»l)»UI  I066S1863 
FUul  I*  LFN  SCCRuf  STATEMENT  -  1FMS!  - 


1*11)  66bSl SbV 

88  »MI  »XMUR2Al)l  MAOS  IRMOI  1  I  •  Ut  II- UA-C  (KPS  1  I  *U  I  MPS  I  I  I 

(.0  10  too 

loci  IF  (A8S1UU)  -UI2I  I.Ll.3.CCi*0»ll  I  CC  -til  1004 
IFtUIMPStl  .61.  20.100  TC  100) 

IT  UR  8*  l  _ 

ITtRA.l 
CO  10  26 
1(  Cl  I T  CM  0*  3 
IT  UR  A*  3 
GO  TC  60 

10C4  KU*WKH  11*0111 

ROM*RHG<MPSl}*MMPSn 
XluT  *X* IR UN  I *RUX  T I /900.*l.t-04 
CO  99  1  * l  .MP3 1 
99  X“0<  I  I  *XMUf 
Go  TO  100 
26  CUM  INUF 

COM*  C. 9*1 KMCl  1I*U(  1I*RHC(MPSI1*UPFSI)  J 
CO  831  11-1. MFSI 
1*  II 
IM  *  I  1-1 

IF(RH0C1I*0(II.GT.CUMI  1  *HPSl— l  1*1 
IMP  HOI  1 1 *UI  1 1 «GI. OOH»  IM=l-l 

IF(RI-Om«UlII-OUrtl  831,832.832  • 

8  It  COM  INCE 

HJ2  RoAL F*YCI l-lYltl-YltMl  I  ♦  IRHC  ( I  I  *1 C 1 1  -DUMl/  (RHO(  t  )*U(  1 1  -RHOI  IM  )«U(  I  ' 

IM  I  / 

GO  TC  733  / 

oC  CoM  INUF 

Ol.’Ml  *3, 8*  lull)  ♦()  IMPS!)  I 
CUV2  *  J  .99*UI  1  1*0.0  1»U  MPSI  I 
CO  721  l»2.MPSI 

IMIUI  1-1  l-l)U*>ll*l  U  1  l-CLMl  1. 1 b  .C.C) 11  =  1 
IF  (I  Ul  l-l  I-DUP2I  *1  U(  1 1 -COM2 )  ,l£  .C.C)  12*1 
7)1  COM  INUF 

80NF  *Y(  11  l-IYIll  I  —  V  C  I1-1))*IU1  1 1 1-CUP1 1  /  1U  HI  T-lTlTr-  IT) 

81  A(.*YI  I2I-IY1I2I-YI  12-1 1 1  AIM  I  2I-CUP2I/IU112I-U1 12-1 1) 

RUAl  F*  AtlSI  8C7.F-R  Tad 

7  I  '  CO  8  31  I*  t.vPSI 

8  11  X'«u*  !)*XMUK2ARhALF»8MCIU*U(l> 

IF  I  JhAll  a  (k  i-AIF- Y(  MPSI  I  I  .01  .C.O  )C  Al  L  EXIT  _  _  231 

X*oT  *XMU(  II 

102  k F  TORN  J60S1992 

£Ni;  366S1603 

FR.)7*  -  b  F  N  9CURCE  STATEMENT  -  IFMS I  - 


3U98  COT  INt  (XPIIC  066S1944 

Cov^CN  R.x.CPCX.P  lU.Mt.FF  .IIPOUT  .F  .OELPi  1  •  CX  .XMUT.XMUL  .OIL  .XMtIK  1.  X 
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TABLE  A-III.  FLOW  FIELD  COMPUTER  PROGRAM  (Continued) 


IMUK2  ,PRNT  .PCNT.UE,  RHCE  .TE.XP3.RFAIF 

CORMCN  I0UT.MPS1  .MHALF.NINIT.1F  I H  S.NPSI < I  RAGf  •  1  PRESS',  ISTA'RT, ! TUkB 
l.' TURA.IJfcT.NS.MOAY.MONTH.PYEAR .IAX2C.1VAR 

f  'KCN  NTNOIE (21 1 • COCFFPI 5  I.AIPHA  (2l",59J.I' 

1.  '  15S1.HI  21  ,  iS  I  .RALPHA  (21<59),PM(59),KU(59I  ,U < 5 9 1 . XLE < 59 1 , T l Tl E (  1 
22..SUM(59).ETA{59» ,XMU(59) «A(59) .AkAlT (591 ,RHCOUTf59) ,U0UT(59 I ,TOU 
3T(  59  ).CP(21.55),HSTG(59),KHCt59),Y(59)  ,XMAX (  3) .  XFRT  (  3)  .CPBARI59 
9lGPA(5S),HBARIS'>),AFITS(2l.2.?l,IF!TSI21.A)  .SPEC  10(21 ) 

(Fit  AX2U.FO.OI  I)X«CCIPS(»0ELPS(  « S (OHAI 11 /XMU (1)  /  XU  ( 1 1  / 12. 0 

IF  ( 1  AX20.EQ. 1 )0X*DEEPSl**2*SIGMA (1 |/XHU{TT7XIE(  IT/WHOI  nVLTTT/A'VC 


CO  1C  I*2,NPS1 

CIVt$*A(i*lUA(I-l)*At  II.Atl) 

DELX«0ELPSI**2*S1GHA(1  l/XLEIll/CIVIS/1.5 
IF  ( 1  AX20.EO.C  1  DEL  X*OE  LX*P S 1  ( ( I 
10  CX*AN!M(UX,OELX) 

"OIFV  i  Y(  2 1  -  V ( 1  >  - 

IF(OX-OIFY)  U2.lt2.163 
163  CX»OlFY/2. 

162  CONTINUE 

OH  ICO  I*2.NPSI 
EX1»CELPSI**2/0X 

"  (F  ( I  AX20.E0.0)  Cxi *£X1  *PSI  ( ( I  . . 

EXlt».5*IA(  IMAIUlll 
FX  12  *. 5*( A(  l)AA(t-ll) 

RUII  l*iFXll*(L(I«l)-U(  1 1 1  +  EX12*  (L  ( l-ll-Ul  1 1 )  )/EXl*U(l  I 

EX4*0.0 

CO  2C  J-l.NS 

20  FX4*EXA*CP(  J»  1 1*1  ALPHA  ( J.1*1)-AIFHA(  J.T-l’fl - 

FX  2*  EX l*CPB AR ( 1 1 

EX 5*XLE( I  )*A(  1  )/SI  GMA(  1 1 

EXft*.5*(EX5  +  XLEt  U1)*A(U1  l/SIGPAIUll  I 

EX7«.5*(EX5yXIE(I-UAA(1-1)/SIGPA(I-11  I 

EX  6*  CP  BAR  (  I  )*A(  I  )/ SIGMA!  II 


C6CS0369 

06050370 


EX 9*.5*(EX8«CPBAR( !♦! I *A(I +  1 I /S IGPAt U 1T1 - 

EX10».5*(fcxe*CPttAKtI-ll*A(l-ll/SlGPA(t-ll) 

EXl3«EXl-EX6-EX7 
EX  14.£X4*FX5/4. 

RT  ( I  ).(U(I*ll-U(t-ll)**2*A(n/EX2/l00151.-23*rtEX9FEXUr>*TrUI»MEX 
1 1J-F  X14  l*T(  I-1UIEX2-E XS-E XIOI *7 1 1 1 1 /  EX2 

on  so  j*i  .n s  - 

40  RALP HA ( J . 1 1 »( E X6*ALPHA ( J , I ♦ 1 )♦ E X X3»ALPHA (.1.1 1 ♦EX7*ALP_HA U,  1-1 1  l/EX 

RT ( l  l«RT( 1 1 ♦CX*JPOX/RHC( 1 1 /CPBAFI 1I/25C37. 807 
RU(II-RU(I)-CX*OPOX/RHO(ll/Um  ’ — — . .  “•  ‘ 


066S2000 


366S20O2 

066S2003 

066S2G0A 

066S2025 

0X6S202B 
G66S2029 
066S20  30 
066S2031 
066S2032 
066S2033 
T55SZ034 
066S2035 
066S2036 
066S2037 


ICO 


CONTINUE 
IFUAX2D.E0.CI 
IF  (IAX2U.E0.il 


066S2060 


FX16*A.0*XMU(ll«C»/CElPSI7DErm- 


200 


E XI 6*2. 0*XMU( l) *CX»RHC ( 1 I *U (1 1/DfclPSI **2 

R'J  ( 1  )*-DX*DPDX/RHO  Il)/U(1)»EX16*  (L~(  2  V-U  t  111  Hjril - - 

00  200  J* l, NS 


RALPHA(J.1I*EX16*XLE(1 1*1  ALPHA! J.2) 
1) 


-AiPHAcj,  n  r/sicpAu  ruownii  i 


RT  ( 1 1 
1T(  1) 


>DX*OP  DX/RHCI  l)/CP0AR(ll/2;C3r.( 


FR02A 


EFN  SCLRCF  STATEMENT  “TF’NTST 


1 


50 


5 


1C 

ICOO 

999 

1001 


CFL*CX  '  " . 

IF ( t  FINIS I5C. l« SC 
IF  IN  I S  *  1 
MINI  T*MPS I 
MHALF*MPSI*MPSI- 1 
CONTINUE 
X»XACX 

CO  3C  1*1.59 

00  5  J*l. NS  *  '  ~ 

ALPHA!  J.I  )*RALPHA(  J,!l 
TI  I)*RTII  I 
U(  1 1  *R  L  ( I  I 

1F(MPSI-MHALFI59S.  1500,1500  . 

IF  ( A 3S I Ut  NP  SI  l-UMPSI  I  l-,001*U(  RFSI I  I  1001 . 100A.  100A 
IFIABSITINPSI  l-TIMPSIl  l-.001*T(FFSII  11002,  1004.1004 


066S2lTf_ 
366S213B 
066S21 39 

066S2146 

<566521 56"" 


C66S2162 

C66S2164 
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TABLE  A- III, 


FLOW  FIELD  COMPUTER  PROGRAM  (Continued) 


u  c;  rxi  i co  j  j*  i.ns 

IMALPHAt  J.PPSII-1  .' -U>  1  C  C  i*  1  CC  J  ,1222  066S2336 

12  22  It  lAHSI  Al  PMAI  J.f.PU  | -ALPHA  (J.HPS1I I- .00 1  *AL  P8A t  J  .HPS 11  1100  3.  1 00  3.  1  _ 


lot  I  CuM  INOI  066S2339 

('«.  10  200o  066S2340 

UC4  “I'M  «HPM»l  G66S234 1 

NP  SI  »MPSI  -  t  066S2342 

C.  TO  2000 

isoo  if  in  i o »o  ■  56'ts irrr 


CALL  FPO/1 
CALI  CCMPut 

IMIftMS.ot.il  HLtuPN 
CElPSl*OClPSI«OEIPSl 
0(1  1600  I  ■  l  .HIM  I 
U(  1 1  *ll(  2»  1-  1) 

Cd  ltSC  J*l.N$ 


ISS 

161 

066S23R8 
066S23B9 
3?>6r2T<>6'  . . 


ln‘J  ALPHA!  J,  |  UALPl-AIJ  ,2*1-11 

loci  ri  *) -r [.>•  i-ii 
“P  >1  MINI  T 
NPSl  *MP;>|-1 
CO  1  70C  I*M|N1T.S1 
-  .  4 1.  J  1*1  IMP  SI  I 
MU  l»l(MPSI  I 
*  00  l  /SC  J  » 1  *N A 

Al  PM  At  J.  I  U  ALPHA  I J  ,MPSi  ) 
U‘0  KAIPHAIJ.  I  I  ‘ALPHA!  J.PPSI ) 
Poll  l«  LIMPS  I  I 
WOO  (Jt  1 1  *<J|MPS1  1 

00  1 80C  I *2 • SS 
1KC0  PSll  IUPS1II  -  1 1*06  LP  SI 
2uL0  Hf  TOWN 


066S2393 

066S2394 

066S2395 

066S2396 

■■'5565Kinr 

066S2403 


066S24G6 

C66S2407 

TJ66524W 

066S2409 

066S241Q 

C66S2411 


toe 


066S2A28 


FMI2S  -  i t.N  M'LRLE  STATEMENT  -  IFMSI 


SUrtP 00  TINE  INI TAL  O66S0651 

L  l«MCN  «.  X.OPCX.POLT.HE  .PE  rOPOU  1  *f  .OSL  PSI  .  CX  .XMUT.XMUL.IIEL  .XNUK  l.X 
1MUK2.PRNT  .PCM.Ut.RHlt  .IE.XP3.kMLf 

Common  lour  .mpsi  .mhalf  .pim  t.i  f  ims.npsi  , i face.  1  FRtss.TsT apt,  itUrb 

1.1  TUPA.  IJE  r  .NS,M0AY.P0MM.PUA8,|AX2C.IVAR 

COMMON  mTMUIE  121).CCEFFP(S»,ALPHA!?l,S9F7MCut09|.RIJCTT1591.RTlS9l  ' 
l  •  T  1 S  91  ,H(  21  .SSI  .P  ALPHA  I?  1.  S9 ).  PS  1 1S9 )  ,RU(  S9)  .01  SSI.XLE1S9I.TITLEI1 
2.21  ,SUM(S91,fTA09l  ,XMU(S9I  .  A I S9 1 .  AkAI  T I S9I  .RHODUT't  STf.UOUTt  597  .  TUU 
IT!  S<»  I. CPI  21  .SS|,HSTGIS9»,RHGIS91  ,7(591  ,XMAX(  3)  . XFRT  1 3 1  ,CPB ARI  59 1, S 
4|t'.PA|59I.HriAK(«9).AF  IT  SI  21.2,7 1  .IF  ITS!  21  .4)  .SPEC  10(211 
XP  »•  C.  C 

ITLRA>|A0S(  ITLRhI 

crucx  066SC680 

IF  LAC  1  *0 

<):>  220  Wl.MPSI 

KII“A  *0.0 


1C 

14 

ICC 

1 1 1 
1 J 

15 

2.0 


22- 

29 


no  ic  j»i.*(S 

i;0“A»DUMA»4lPHA|  J,  II 
II  IDUMA.EG.O.C1UO  TO  2  2C 
IF IAPS(OUMA-t .  I- .01 1  1  1,  13.  ICC 
WP  m  I  6.1  11  ICLPA.I 
It  LAC  l  *  1 

Fo-MAt  I  13.(0  SIGMA  ALPHA  >lPtlS.7,  12HP0INT  NOPHEK  I  S  ,6MlN  PSll 

CO  IS  J.l.NS  . .  . . . 

Al  ?mA(  J.I  IMIPhAI  J.l  l/OLPA 
COM  ISLE 

If  I  I  FI  AO  1  .10.  I  I  CALL  I  X|  I 
(X)  228  1*2.  SS 
XI  FI  l)*XI  U  II 

SIO'AI  I  I  =  S I  r,M  A  |  I  1  . 

PS  II  1 1  *PS 1 1  I-  1  ItOtLPSI 

nu“MY=3.ost-g 
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__  066SC69_6 _ 

O66S0698 

C66SC699 

29 

066S0700 

U6S0701 

l;66S'07<5? - 

066S0703 

066S0704 


P*  Ci. tr  FPI  1 1  ♦XMLCiFFPI  2UX*IC(.tttri  JI*X»|CCtFFPI4)tX»CCEFFPI5!  1 1 1 


TABLE  A— III.  FLOW  FIELD  COMPUTER  PROGRAM  (Continued) 


on  vc  i«mpsi,ss 

RT  ( I  )*  T  IMPS I ) 

- rm»T(Mpsn  ■  " 

00  SO  J-l.NS 

RALPHAU.II-ALPHAt  J.MPSI) 

80  ALPHA!  J.ll-ALPHACJ.MPST) 

RUtl  )*U(MPS l ) 

SO  U(  1 1 -UIMPSl ) 

*“■  «RItE(4'.3An<TITLEU>.l*Ul2»‘  ” 

IF  ( I  AX2t).  EQ  .C  1  WRI  Tb  (6.616) 
lfllAX20.EQ.il  WRI TEI6.61?) 

HR  IT E( 6.888  I  RCNTH.MDAY.MYEAR 
WR  II  £1 6,999  I 
Aft  IT  El  6.609) 

- -WUTWiMtr  •  •  . . 

WR  IT  El  6,6 12 1 
HR  IT  E(  6.6 15  I 
UR  IT  El  6.6 18 1 

IF  1 1  PRESS-2 1  tCl.6C2.6CA 
601  URITEI6.7C1 1CCEFFP 1 1 1 

-  GO'TOiCA  •  . . 

6C2  HR  IT  El  6.702  I (CCEFF  PI tl .1*1.6) 

FR02S  -  EFN  SCURCE  STATEMENT 


066S075A 

066SC755 

”’556SoT?l 

066S0758 

C66S0759 

066S0762 

066S0766 


5F?CR!T5~' 


-  IFMST  “ 


ECCXONT'INUE  . 

WHITEI6.999) 

HR  IT  El  6.806 1  SIGMA!  It.  XL  Ell) 

HR  IT  6(6,969  I 

HR  IT E (  6,809 )  ITliRA 

UR  IT  FI  6,8 10  IXMLK 1,  XMUK2 

- XHLK  1-XMUK1*DUMMY  0ST6SC66E 

HR  1TEI  6,999) 

00  171  1-1,3  " 

171  HR  IT  El  6.812  I  XPR  Tl  1 1  .  XMAX  1 1  I 
HR  IT  El  6,999  ) 

HR  IT  El  6,6 1.3  )  UIMPSI )  ,T  IMPS! )  . 

*  RETURN  ^  - -  ——————  — 

3A1  FORM AT  1 1H1 , 2AX  .3  3HGENE  RALI 2EC  FRC2EN  MIXING  PROGRAM/ 1H0.5X.12A6I 
8E8  FORMAT  1 1H0,  3AX  .AHOATEI 3. 1H  /  .  12 «  ip/  ,12 1  '  " 

6C9  "FORMAT 1 18)1*  VbLOC  ITY—IFT/SEC)  ,20X .'17HENTHAI  PY-IBTU/LBI) 

611  FORMAT  1 29H  TEMPERATURE -(DEGREES  KELVIN)  »9X»'2ZH|T$'I=fC$l0C-S7$6O**l/' 

121  )  _ 

SH!  To  RM  AT  1 22H  CEASITY-ISl  LGS/FT**3 1  ,ifX.l^W(rE5SURF-aH/n,«21T 
613  FORMAT!  19)1  EDGE  VE  LOCI  TYIUE  l*E 1 1 • A »15X .21HE0GE  T  EMPERATURE I TE ) « El  l 

l.AI  . 

615  FORMAT  1 29H  VI  SCO  SI  TY(K0)-ILB*SEC/FT**2  )  ,9X,12HRACI  I-IFEET) ) 

6T6  FORMAT!  33X.17HAXI SYMMETRIC  FLOH)  . .  "  ' 

617  FORMAT!  32X.2CFTHO  DIMENSIONAL  FICU) 

618  FORM  ATI 25H  SPECIES-IMASS  FRACTICNSJ) .  . 

8C8  FORMAT! 16H  PRANDTL  NUMUER-F5.2 , 2AX .13HLEUIS  NUMBER-F5. 21 

809  FORMAT!  18H  VISCOSITY  OPTION-12)  . . 

810  FQRMATI31H  LAMINAR  VISCOSITY  COEFF IC IENT-F5 . 2,9X  .32HT0RBUL  ENT  VISC 
10SITY  COEFF ICIENT-F7.A) 

812  FORMAT  1 12H  PRINT  EVERYF7.3.5H  FEE  T  ,2 IX , SHUNT  1 L  X-F8.3.5H  FEET) 

ST99 'FORMAT!  1H0)  '  '  . 

701  FORM  AT  1 1H0,  20X,23HCONSTANT  PRESSLRE  P-E15.7) 

702  FORMAT!  1H0.  3CX.  63HPULYN0MI  AL  PRESSLRE  FIT'  >»TaB'*IXT*C'*TX**2TfD«066SC960 

1(X**3I*E*(X**A)/15X,2HA»IPE15.7,3H  B*1PE15,7,3H  C-1PE15.7.3H  0- IPE066SC961 
215.7, )F  E-1PE15.7)  066S0962 

END  066S1008 

FR026  -  EFN  SOURCE  STATEMENT  -  IFMS)  - 


SUBROUTINE  CQMPUT 

DIMENSION  H  JACK!  60,251  ,P  JACK  160 ,25 1 ,  THAT  EM  (60, 25)  »XJ"ACK(60  )  ' 

CO  MM  CN  R, X, OPCX, POUT ,HE, RE  «OPUU! , P  tUELPSI ,CX ,XMUT,XMUL,0EL .XMUKl.X 
IMUK2 ,PRNT .PCNT.UE.RHOE ,TE,XP3.RFALF 

COMMON  IOUT.MPSI ,MHALF .MINT  T,I F  I M  S.NPS I , I  PAGE, I  PRESS, 1ST  ART,  ITURB 
I.ITURA,  UET , NS, MDAY, MONTH,  MYEAR.I  AX20,'IVAR 


53 


TABLE  A— III.  FLOW  FIELD  COMPUTER  PROGRAM  (Continued) 


COMMON  MMULEIil  I  .COEFFPI5  >. ALPHA! 21 ,59)  ,H0UTI  59  >,  ROOT T l  59  ) , RT  (  59) 

1,1(591  ,H(  21 « 55)  tKALPHA  ! 2 1 <  59 ) , P S I  (59 ) , RU( 59)  ,U('59)'VXCE’( 59)~TtTL'£n - 

22)  .SUM!  591,  tT A!  591  .CMUISO)  ,  A ! 59 )  ,AKA  1 T 1 59 )  , RH00UK59)  ,UQUT  159) ,  TOO 

3TI  59  ).CP(  21,5SI,HSTG(59),RHC(59>,Y(5S)  , XKAX  1 3)1  XPRH 3 )  ,'CPBW!  5'9  )TS - 

5IGPA(5S),H0AKt59l,AF|rs<21  ,2,7)  , IF  1  7 S ( 21 , A ) , SPEC  10121) 

loot « tour  ♦!  - - 06550075 - 

IHXP3.E0.0.CI  XP3  =  XPR  Till 

IP  ( I  AX  20,  CO  •  C )  GO  ro  52  -  --  .  .  -- 

00  51  1  =  1  ,MPS  I  « 

51  At t ) =XMU( I ) *RHC( I ) *U 1 1  )  .  . 

GO  TO  53 

52  A!  II  *0.0  - - - 

A(  21  =OELPSI  *XML(  21 

00  55  1  *3 , MPS I  - - - -  -  '  ’ 

55  A (  !)*XPO(  I )«RHC( 1) *U(  I )  *Y( 1 1  *Y  I  I  l/PSltt) 

53  CONTINLE  - - - -  ... 

RflAk  CL*X 

IP  ( t  FINIS, EC,  C)  GO  TO  50C  - 7 - - - 

CO  500  1=  1,3 

IF  (R  8AR0L  ,L  T .  XNAXl  1 1  )G0  TO  501  '  . .  ' 

500  CONTINUE 

GO  TO  888  -  . .  . . 

5C1  PRNT  *XPRT  (  I  I 

IF  (R8AR0L  ,LT.  XP3 )  RETURN  ‘  - 

00  6695  1*1,3 

IF  (R HAKOL  , L  T ,  XHA X (  I ) )  GC  TC  668e  . . 

6659  COAT  INLE 

888  IF  IN  IS *2  •  ’  - - - -  ‘ 

GO  TO  555 

6639  XP  3=  XP  3*XPR  T(  I )  - 

6C0  CONTINUE 

5  55  COM  INUE  . . .  . . 

572  F0RMAT(E10.8, 15) 

IFIX.EC.O.)  R  EAD(  5 » 572  )  ROC.AORLN  - - -  •'  '  ~68 

IF  (X  ,E0  .0  ,  )  I  XP*C 

IFIX.EC.O.)  GC  TO  eC2  - - 

CHGX  *X-PREVX 

CM  PR  =XPRT  1 1  )-C ,  1  - — . . 

IF(CHGX.LT.CMPR)  RETURN 

8C2  PR  EV X*X  '  - - - * -  ~  *s —  — - - 

IXP-  IXPM 

IFIIXP.E0.61)  GO  TC  999  - - “■ 

TT  ES  T*  1 .0  5*  I  ( PPS I  ) 

XJACKI  IXP  ) * X  - - - - * . . 

CO  525  KUNK*1 ,26 

RJACKI  [XP  , KUNK  )>  Y(  KUNK )  - * - - 

SJ  ACK*C.O 

CO  530  KONK-l.NS  - 

530  SJ  ACK»SJACK ♦ ALPHA! KONK  , KUNK I /HTPCIE I KUNK ) 

FR0Z6  -  EFN  SOURCE  STATEMENT'"  -  IFTT5T  - 


-  PJACK!  IXP  ,K  UNK  )*  AL  PHA(  3  .KUNK)  /  W  TPtlETi  fTJXTEK  : 

TR  AT  EM  I  IXP,  KUNK  )*T  (KUNK  I 

IFtTRATEM!  IXP, KUNK). LE.TTEST)  GC  TC  550" 

52  5  CONTINUE 

550  IF  (T  (KUNK  I.GT.TIMPSI  )l  GU  TC  555  <  , 

555  PJACK!  IXP, KUNK  )=C.C0033  9 

OYP.Y(KUNK-l)  ♦  Yl  2)  ‘  - 

GO  TO  558 

555  kTRI  »ALOG  (  T  (KUNK- 1 I  /  T(  MPS1  I  -  l.)  - 

n'TR2  *ALUG  ( T  ( KUNK  1/ T ( MPS!  I  -1.) 

WTRP  *ALOG (0,05)  . . . 

YT  R1 *YI KUNK-1 ) 

. YTR2  *Y  IKUNKl"  '  - - 

YTRP  *(  YTP  I«YTRI*(WTRP-V,TR2  J5YTR25YTR2* (HTR 1-RTRP ) )/ 1 WTR1-HTR2  ) 
OY  P=  SORT ( YTRP )  - - 


2  T  Rl *AL0G I P  JACK! IX  P , KUNK-1 )— C.0CC33)  128 

2TR2  *ALUG IPJACKI IXP, KUNK )-C,C0033  )  136 

BTRA  *(  2TRI-2TR2)/!  YTR2*YTR2-YTR  1 «  YTR1) 

CTRA  =2  TR1 *B  TRA*YTR 1* YTR1  . . 
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TABLE  A— III.  FLOW  FIELD  COMPUTER  PROGRAM  (Continued) 


PJACK1  lXP,KUNK)s0.0003  3*r.X?(CTHA-6TRA*0YP<-CYP) 

568  RJACKI  IXP  ,KUNK  )*DYP 

TP  AT  EM  ( IXP , KUNK)  *T  Tk  ST 
lFlKUNK.fcQ.25)  GO  TO  565 
KUNK PI aKUNK ♦ 1 

00  666  1JK=KUNKP1, 25  ~  '  *  ' 

PJ  ACM  IXP : 1 JK )  *P JACK!  I  XP.KUNK) 

RJACKI  IXP.tJKlsKJACKUXP.KUNK) 

566  TRATEMUXP,  I  JK  I  =T  TE  ST 

565  IM1XP.N0.60)  GO  TC  56  7 
ATR=0. 

CTRsO.  '* 

00  566  UK.s2C.63 
TY  sR  JACK!  JK,  251-RGO 
ATRs  ATR*T Y 

566  CTRsCTRKX JACK  1 1JK) 

CTRsATK/CTR 

DO' 568  1JK*1,6C  ! 

568  RJACKI  1  JK  ,  25  )  sROC+OTR*X  JACK  (  I  JK  t 
00  573  NP « 1 »6C 
DO  573  iP*l«26 

If  (TKATfcMINP,  IPI.EQ.TTKST)  R JAC K  l  AP, 1 P )  =  R J ACK 1 NPY25 >  *  ~ 

5  73  CUNT  INUE 

00  5.76  NP*2 « 6C 

If (R JACK1NP  ,25) , LE .R JACK (NP-1 » 2 5 ) )  R JACK  I NP, 25 ) s P JACK! NP-l ,25 ) *0. 1 
5  76  CONTINUE  “ 

HR  IT  El  6,9000 ) 

90  CO  FORK  ATI  1H1)  . 

90 Cl  FORMAT  1  1H  ,6E12.5I 

00*  5  76  NP*1,6C  “  : 

5  76  HR 'T El  6,9091)  XJACK1MP  )  ,KJ  ACK  INF  ,25)  ,  TRATEK( NP ,2  5 ) ,  PJAC.K (NP,  2 5 ) 

11=25  . 

00  577  NP* 1 , 6C 

XJACK1NP)*XJACK1NPI*30.6E  ‘  . .  ' 

00  6 7H  KK*1 ,25 

'5ra  RjACK(NP,KK).-.RJACKINP,KK)*30.6fi 
5  77  CONTINUE 

FR026  -  FFN  SOURCE  STATEMENT"  -  “IF'MTT  - 

*  "Cl  MENS  ION  CAT(25),RAT(25),PAT(25>  “ 

00  750  NXP*  1, 6C 

00  751  NRF-*  1, 25  .  . .  . . "* . * . 

lflTRATcMNXP.NRPl.EQ.  TTEST)  GO  TC  752 

751  CONTINUE  ‘ 

752  RAPAXsRJACKINXP.NRP) 

'NATPLsNRP 

DO  753  NN»1 ,NRP 

IF1RJACK1NXP.NN)  .GT.RAPAX)  GO  TC  756  * 

7  53  CONTINUE 

GO  TO  755  " 

7  56  NRP*NN 

RX*I  RJACK I  NXF.NRP-  1)*R  JACK  (NXP.NRPAlll/Z,”" . . 

PTl*(TRATEM(NXP,NRP-l)-T(MPSI) )  / 1  TRATEPI NXF.NRPa'I )-T  1 MPS I ) ) 
IF1PT1.LT. l.C)  PT1«1.0 
PT  2*~ ALCG 1 P  T 1 ) 

PT  3*RX*RX-R  JACK! NX  P, NRP- 1 ) **2 

PT 6* RJACK  1  UXP.NRP6  1 1  ** 2-RJ ACKl  NXF  ,NRP-1 )  **2  __ 

PT5*  1PT3/PT6)*PT2  *  "  .  “  * . 

PT  6*  TRAT£M(NXP,NRP-l)-T(PPSl) 

TR AT  EMINXP.NRP  )*  T1  PPM  )6PT6*DXP  1PT5I 

PT  1*  ( (  PJAC.K  1  NXF, NR P- 1 )  -C .00033 )  1/  (PJACKINXF.NRP*  ll-U.00033) 

IF  (PTl.LT .1.0)  PT1  *1,3 
PT2S-AL0G1PT1) 

PT5=  1PT3/PT6)  *PT  2  ' 

PT  6*  PJAUK I NXP.NRP- 1 1-0. CC033 

PJ  AC  K( NXP ,NRP ) *0 ,000336PT6*EXP 1 FT5 ) 

RJACKt  NXP, NPP)«RX 

755  CONTINUE  . 

OELTRsR AM AX/26. 

CO  756  NRP«2,25  *  ' 


J. 000331 
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TABLE  A-III.  FLOW  FIELD  COMPUTES  PROGRAM  (Continued) 


CNT*NKP-1 

RA  71  NRP  )*CN  T*DEL  TR 
PI  -O.OOOl  I 
Tt  MIMPSI  ) 

CO  761  IJK*2.NATPL 

IF1RAT(NRP).LT..KJACK<NXP,IJK)I  GC  70  75!T 
7  67  CUNT  IN  UP' 

7  56  R 1  =R  JACK!  NXP » I  JK-l ) 

R2=«  JACKINXP.  IJM 
RX=RATCNRP) 

Tl*TKATtM<NXP.IJK-t) 

T2=TRATFM<NXP.IJK) 

PI  =P JACKINXP.  1  JK-l  I 
P2=P  JACKINXP,  1JKI 
PT  |=<KX*KX-R1*RI  1/  <R2*R2-Rl*Ril 
PT2=  C  T  t-T  1»  /IT2-TI  I 
IFIPI2.LT. 1.0)  P  72  »l .  C 
PT2=-AL0GIPT2)»PT1 
CAMNRPI*  II  ♦171-71  )*EXP(PT21 
PT 2* IP  l-P  11  /IP2-P!  ) 

IFIPT2.1T.  l.C)  P  72*1,0 
PT  2*-ALOG(PT2)*PTl 
PAT(NRPI  =  PU(  PI-PI  )»EXP(PT2) 

7  56  CONTINUE 

0.1  759  JK=2,2E 

FPIlZtj  -  £FN  SOUHCF  STATEMENT 


RJACMNXP»JK)*RATI  JKI  .  ' 

PJ  ACM  NXP  «  JK  I  *PAT(  JK| 

759  TRAT  tMlNXP.JK  1  *C  AT  (JK) 

WRITE!  IOMTIHF.I  It  ,1*1. 121 
WRITE!  10)  flljRLN 
GO  579  .NP *  1  •  6 C 
WRITE!  10)  XJACKlNPI.il 

wR IT  El  10) (TkATEMINP.JJI .RJACKINF.JJI , P JACK ( NP, J J ) , J J* 1 , 
5  79  COMINLE  '  ~ 

WR  I T  El fc « 766  )  XJACKINXP) 

DO  760  JK*1,25 

7  60  wR  ITEI  6.707  IR  JACK!  NXP,  JK  I ,  TP  ATEH  NXP,  JK)  ,P  JACK  I  NXP,  JK) 

7 50  COMINLE  ‘  ~ 

7E6  FORMAT  (  1H  ,E1!.A) 

767  FORMAT!  1H  ,2615. A» 

567  l)J  75  11*1, MPSI 
I*PPSI  +  1-I1 
70UI II  I  *T  I  I  )/l  IMPS  I ) 

UOLTII  )*UII  l/UtPPSI) 

RHC0L7  I  I )  *KHOI  I )  /RHOI.MPSI I 

HS7GII  )*HHAR(  I  )*-U(  I )  *0 «  D/50C75.61A  "  . 

FOOT  1 1  )«HSTG< 11/H$T0(MPSI) 

SUP(  11*0.0 
CO  1  l  J*).,NS 

tl  SUPI  II*S0flC  ll+ALPHAU.!)  *  ' 

25  CONTINUE 

IPAGE* IPAGE+l 

WRITE!  6,201)1  I1TLEII  I .  I  *  1 , 12 1 . 1  PAGE 
WRITE(6,721)  NCR  UN 
721  FORMAT!  1H0.  1CFRUN  NUPHER.l  A) 

WR  IT  El  6, 102 )  X.UX.IGUT  “  ‘  " 

WR IT  El  6,600 ) XPLL ,XPU I , RHALF 
WR  ITEI6.609)  F.OPOX 
wR  1 1 E 1 6 , 1 C  7 1 

C.l  10  1*1. MPSI 

10  WR  IIFI 6,207)  I  «LJUT  1 1 ),TCLTII),RFCCLT|1) ,HOUT ( I ) , SUM! I ) , 

‘  '  1.1  '  - 

K=M|  NO!  7,  .NS) 

WR!TF(6.20ll(TITLE(I),l«l,12l,IFACf  -  -- 

wR 1 1 E ( 6 • 1  OH )  (SPEC (Oil ). 1*1, Kl 

00  20  1*1, MPSI 

20  WR 1 7  E ( 6 , 20H I  I , I AL PHA I J , 1 1 , J*1 , X ) 

IFlNS.Lt.7)  RETURN  . . . 


'06«C0fr~ 


066S0056 

‘KZitasr 

066S0058 


ort'Scl'aT' 
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TABLE  A-III.  FLOW  FIELD  COMPUTER  PROGRAM  (Concluded) 


K-KIMMIA.NSI 

hRiTe(6.2onnmem.i«uizi,iFAct 
WRITE!  6. 108)  (SPEC  10(1  >,1*8,1'.) 

00  3A  1*1, MPS! 

35  HR  IT  El  6,208  )  l , I AL PHA t J , I > , J*H , K > 

lRN'S.'CE.TA  )  "RETURN  '  ■"  ~ . . 

K* Ml  NO!  21  ,NS ) 

‘  wRtTEl6,201)UiritM).I*l,12),lFAGE 
VtR  IT  El  6 *  1  OR )  (SPEC  10(1  ), 1*15, K) 

CO  "35  l*l*MPSt 

35  NR  IT  El  6,208  I  1 . 1  ALPHA!  J, 1 1  .  J*15  ,K ) 

— vnufiH —  -  -  - - - 

999  IF  IN  IS«2 

"  FRUZ6  -  EFN  SOURCE  STATEMENT" 


A  3C 
A  3  7 

AA  7 


A5t 

A63 

"  i 

a73 


-“irKCSi — 


- Wen/ftST  '  * . — - 

2GB  fORMAT 1 13, 7E1 1,3 ) 

108  FORMAT(3HOPT.7(3X,A6,2XI) 

2C1  FORMAT! 1H1, 3X , 12A6.2X) AHPA0EI3 ) 

207  FORMAT (13«?£15<6<1A)  ‘ 

10  7  FORMAT  I  3H0PT  «6X»  AHU/UE  ,  11X  ,AHT  /  TE  ,9X  ,8HRH0/RH0£,9X*  AHH/HE,  1QX »  5HS I 

- 1T.M AVI 2 X.lHY , 1 3X, 3HPS 1 , 7XV2HPTI - 1  " 

608  FORMAT! 6H0MU  L*E  12 . 5 , 7  X, 5HMC  T* E 1 3 . 5 , 1 3X ,6HRhALF*E 1 3 .5 ) 


~6C9“faHMAn  3H0P «E  11.  3, 10X,  5hCPCXiE  1 1  .'3  T 
102  FORMAT! 3H0X*E13.5, 10 X, 8H0ELTA  X*E  13. 5, 10X , 6HSTEPS*! A > 
.  ENO  ' 


"3W5WTB 
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